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We introduce two probabilistic models for N interacting Brown- 
ian motions moving in a trap in M. d under mutually repellent forces. 
The two models are denned in terms of transformed path measures 
on finite time intervals under a trap Hamiltonian and two respective 
pair-interaction Hamiltonians. The first pair interaction exhibits a 
particle repellency, while the second one imposes a path repellency. 

We analyze both models in the limit of diverging time with fixed 
number N of Brownian motions. In particular, we prove large devi- 
ations principles for the normalized occupation measures. The min- 
imizers of the rate functions are related to a certain associated op- 
erator, the Hamilton operator for a system of ./V interacting trapped 
particles. More precisely, in the particle-repellency model, the mini- 
mizer is its ground state, and in the path-repellency model, the min- 
imizers are its ground product-states. In the case of path-repellency, 
we also discuss the case of a Dirac-type interaction, which is rigor- 
ously defined in terms of Brownian intersection local times. We prove 
a large-deviation result for a discrete variant of the model. 

This study is a contribution to the search for a mathematical for- 
mulation of the quantum system of N trapped interacting bosons 
as a model for Bose-Einstein condensation, motivated by the suc- 
cess of the famous 1995 experiments. Recently, Lieb et al. described 
the large- A behavior of the ground state in terms of the well-known 
Gross-Pitaevskii formula, involving the scattering length of the pair 
potential. We prove that the large- A behavior of the ground product- 
states is also described by the Gross-Pitaevskii formula, however, 
with the scattering length of the pair potential replaced by its inte- 
gral. 
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1. Introduction and results. 

1.1. Introduction and motivation. Consider a large number of identical 
quantum particles in a trap under the influence of a mutually repellent pair 
interaction. We assume that their wave function is invariant under permu- 
tations of the single-particle variables. In physics particle ensembles whose 
many-body wave functions have this invariance property are called boson 
systems. A specific property of bosons is that, for large number of parti- 
cles, at extremely low temperature, they undergo a phase transition, the so 
called Bose-Einstein condensation. This formally means that a macroscopic 
portion, the condensate, is described by one single-particle wave function. 
Bose-Einstein condensation was theoretically predicted by S. N. Bose and 
A. Einstein in the 1920s, and a huge theoretical literature has been accumu- 
lated since. After appropriate cooling methods had been developed, the first 
experimental realization of this condensation succeeded in 1995 [4, 6, 10]. For 
this remarkable achievement, the Nobel prize in physics 2001 was awarded 
to E. A. Cornell, W. Ketterle and C. E. Wieman. A comprehensive account 
on Bose-Einstein condensation is the recent monograph [28]. 

Motivated by the experimental success, in a series of papers Lieb, Seiringer 
and Yngvason [22, 23, 24, 25] obtained a mathematical foundation of Bose- 
Einstein condensation at zero temperature. The mathematical formulation 
of the iV-particle boson system is in terms of an A-particle Hamilton oper- 
ator, H.n, whose ground states describe the bosons under the influence of a 
trap potential and a pair potential, see (1.7). Lieb et al. rigorously proved 
that the ground state energy per particle of "Hn (after proper rescaling of 
the pair potential) converges toward the energy of the well-known Gross- 
Pitaevskii functional, and the ground state is approximated by the TV-fold 
product of the Gross-Pitaevskii minimizer. Moreover, they also showed the 
convergence of the reduced density matrix, which implies the Bose-Einstein 
condensation. As had been generally predicted, the scattering length of the 
pair interaction potential plays a key role in this description. 

Much thermodynamic information about the boson system is contained 
in the traces of the Boltzmann factor e~^ N for j3 > 0, like the free energy, 
or the pressure. Since the 1960s, interacting Brownian motions are gener- 
ally used for probabilistic representations for these traces. The parameter 
/3, which is interpreted as the inverse temperature of the system, is then the 
length of the time interval of the Brownian motions. However, the traces do 
not contain much information about the ground state. Since the pioneering 
work of Donsker and Varadhan in the early 1970s, it is basically known that 
the ground states are intimately linked with the Brownian occupation mea- 
sures. This link is established via the theory of large deviations for diverging 
time, which corresponds to vanishing temperature. 
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In the present article we study models of a fixed number N of trapped 
interacting Brownian motions given in terms of transformed measures for 
paths of length j3. A trap Hamiltonian, putting "hard" or "soft" walls, keeps 
the motions in a bounded region. We consider basically two different types 
of pair-interaction Hamiltonians, both imposing mutually repellent inter- 
action: the first one, which we call the canonical ensemble model, imposes 
pariic/e-repellency, while the second, which we call the Hartree model, im- 
poses pai/i-repellency. More precisely, in the canonical ensemble model the N 
motions interact with each other at common time units, and in the Hartree 
model, the paths interact with the mean of the paths of the other motions. 
The canonical ensemble model comprises the above mentioned standard rep- 
resentation of the traces of e~^ N . The Hartree model is apparently intro- 
duced here for the first time; it cannot be represented in terms of the trace 
of any Boltzmann factor. 

It is the main aim of this paper to study the large-/? behavior of the canon- 
ical ensemble and the Hartree model in terms of large deviations principles 
for the Brownian occupation measures. One of our main results is a large 
deviation principle for the joint occupation measure of the iV-tuple of the 
motions in the canonical ensemble model (Theorem 1.5). The rate function 
is given by the energy of the operator 7~Cn- In particular, we prove a law 
of large numbers for the occupation measure toward the ground state. The 
large-/? behavior of the Hartree model is different and is described in terms of 
the ground product states of 7Y/v, that is, the minimizers of the energy of Wat 
among all product states. Our main result here is a large deviation principle 
for the iV-tuple of the occupation measures of the motions (Theorem 1.7). 
The rate function is the T^Ar-energy of the product of the components of 
the tuple. In particular, we have a law of large numbers for the occupation 
measure tuple toward the set of the ground product states. We also discuss a 
mathematical idealization of the Hartree model, where the pair-interaction 
potential is replaced by the Dirac-measure at zero. This model (which is 
trivial for the canonical ensemble model) is defined in terms of Brownian 
intersection local times, an object whose large deviations properties are cur- 
rently much studied from a probabilistic point of view. We prove analogous 
large deviation results for a discrete variant of the model (Theorem 1.12). 
Finally, we also study the ground product states of Wat in the limit N — > oo 
(Theorem 1.14). In dimension d = 3, the rescaling of the pair potential is the 
same as was considered by Lieb et al. It turns out that the energy converges 
toward the Gross-Pitaevskii formula, and the mean converges toward the 
minimizer. Instead of the scattering length of the pair-interaction potential, 
its integral, which is a strictly larger number, is the decisive parameter. 
In dimension d = 2, we prove the same result, however, the scaling is here 
different from the one considered by Lieb et al. 
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The replacement of the ground state energy by the ground product state 
energy is known as the Hartree-Fock approach [12]. One of the main novel- 
ties of the present article is a probabilistic version of the Hartree-Fock idea 
for interacting Brownian motions; indeed, the product states are interpreted 
in terms of the occupation measures in the Hartree model. Another main 
novelty is a rigorous proof of the asymptotic relation, for diverging particle 
number, between the ground product states and the Gross-Pitaevskii mini- 
mizer with parameter equal to the integral of the pair-interaction potential. 
This relation and the one between (unrestricted) ground states and the scat- 
tering length is a general issue and has been phenomenologically discussed 
elsewhere [14, 29]. However, the former has previously not been rigorously 
proved for the iV-body problem. 

The present paper provides the mathematical basis for the rigorous proba- 
bilistic analysis of interacting Brownian motions as models for large quantum 
particle systems at positive temperature. In [1] we give an asymptotic anal- 
ysis of the Hartree model at fixed finite time for diverging particle number. 
Future work will be devoted to the same question for the canonical ensem- 
ble (i.e., for the trace of e _/3Wjv ), which is important from the physical point 
of view, but also more ambitious. Specific permutation symmetric systems 
of Brownian bridges, characteristic for boson systems, will be analyzed in 
[2]. The latter hopefully leads to a probabilistic model for Bose-Einstein 
condensation at fixed positive temperature. 

The remainder of this section is structured as follows. We first define the 
two basic variational problems in Section 1.2. In Section 1.3 we define the 
canonical ensemble model and the Hartree model. Our results on the large 
deviations behaviors of these two models are also presented in Section 1.3. In 
Section 1.4 we discuss the Dirac-type interaction. The large- N behavior of 
our key variational formulas is treated in Section 1.5. The notion of a large- 
deviation principle and of the scattering length are recalled in Section 1.6. 

The subsequent sections are devoted to the proofs: in Sections 2, 3 and 
5 we prove our results on the Hartree model, its Dirac-type variant and 
the canonical ensemble model, respectively, and in Section 4 we analyze 
the ground product states and its large- iV behavior. Throughout the paper 
we consider only dimensions d > 2; in Section 1.5 we restrict to the case 
de {2,3}. 

1.2. The variational problems. We introduce two fundamental variational 
problems for the energy of N particles in M. d under the influence of two po- 
tentials. 

1.2.1. The potentials. Our two fundamental ingredients are a trap poten- 
tial, W, and a pair-interaction potential, v. Our assumptions on W are the 
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following: 



W : R d -> [0, oo] is measurable and locally integrable on {W < 00} 



(1.1) 



with lim inf W(x) = 00. 
R^oo\x\>R 



In order to avoid trivialities, we assume that {W < 00} is either equal to M. d 
or is a bounded connected open set containing the origin. 

Our assumptions on v are the following. By B r (x) we denote the open 
ball with radius r around x G M. d : 

(12) V '' ^' ~~ > ^ ^ {+°° } i s measurable and bounded from below, 

a := sup{r > : v(r) = 00} G [0, 00), u l[r?,oo) 1S bounded V77 > a. 

Note that we also admit v(a) = +00. We are mainly interested in the case 
where v has a singularity, that is, either a > 0, or a = and Homo v(r) = 00. 
Examples include also super-stable potentials and potentials of Lennard- 
Jones type [30]. According to integrability properties near the origin, we 
distinguish two different classes as follows. 

Definition 1.1. We call the interaction potential v a soft-core potential 
if a = and f Bl r \ v(\x\) dx < +00. Otherwise [i.e., if a > 0, or if a = and 
Jb 1 (o) v(\x\)dx = +00], we call the interaction potential a hard-core potential. 

We shall need the following <iiV-dimensional versions of the trap and the 
interaction potential: 



where x = (x±, . . . ,xn) G Mr . Our potential 2U + D is locally integrable 
precisely on the set 



N 



(1.3) 



W(x) = ^2W(xi) and t>(x) = J] v(\ 





(1.4) 




if v is soft-core, 
if v is hard-core 



where 



(1.5) 



U v := p| {xe R dN : \x { 



Xj\ >??} 



77 > 0. 



l<i<j<N 



Note that is a connected set in d > 2 
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1.2.2. The ground state energy. Our first fundamental object is the ground- 
state energy per particle 

Xn = 4; inf (h,HNh) 
q g\ N hem(n):\\h\\ 2 =V 1 

of the iV-particle Hamilton operator 

(1.7) Wtv = -A + 2U + D on L 2 (Q). 

Here iZ 1 ^) = {/ € L 2 {VL):Vf £ L 2 (0)} is the usual Sobolev space, and V 
is the distributional gradient. We summarize some facts about the ground 
states of TCn- 

Lemma 1.2 (Ground states of Hn)- Fix N G N. 

(i) There is a unique minimizer h* E i? 1 (0) on the right-hand side of 
(1.6). 

(ii) TTte minimizer h* satisfies the variational equation 

(1.8) A/i* = 2U/i* + 0/i* -iVxA^*- 

(iii) XTie minimizer h* is positive everywhere on £1 and continuously dif- 
ferentiable, and its first partial derivatives are a- Holder continuous for any 
a < 1. 

Proof. The proof is standard and is therefore omitted. See the proof 
of Lemma 1.3 below or Appendix A of [22] for the treatment of similar 
problems. □ 

1.2.3. The ground product state energy. Introduce the ground product 
state energy oiTCjy, that is, 

(1.9) x { n ) = T7 inf (hi®---®h N ,'HNhi<8f-<8>hN)- 

N ft 1 ,..,/ 1Jv eiJl(R'i) : 11/14112=1 Vi 

The replacement of the ground state energy, xn, by the ground product 
state energy, Xn > 1S known as the Hartree-Fock approach [12]. Sometimes, 
the formula in (1.9) is called the Hartree formula. Obviously, 

(1.10) X^>XN. 
We can also write 

i r N 

N N hi,...,h N em(®i)-. ||Ma=i vi ■ 



i=l 



+ (hlVhf, 

l<i<j<N 
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where V denotes the integral operator with kernel v o \ ■ \ , either defined 
for functions by Vf(x) = f^dv(\x — y\)f(y)dy or for measures by Vfi(x) = 
J Rd H(dy)v(\x -y\). 

The main assertions on the formula in (1.9) and its minimizers are sum- 
marized as follows. 

Lemma 1.3 (Product ground states oIHn)- Fix N £ N. 

(i) There exists at least one minimizer (hi , . . . , hw) of the right-hand 
side of (1.11). The set of minimizers is compact and invariant under per- 
mutation of the functions hi, . . . , hp?. 

(ii) Any minimizer (hi, . . . , hjy) satisfies the system of differential equa- 
tions 

(1.12) Ah i = -X i h i + Wh i + h i J2Vh 2 j , i = l,...,N, 

with Xi = \\Vhi\\l + (W, hf) + Ej^i( h i> Vh ^- Furthermore, \\hi (U < C d {\i - 
(N -\)\vS.v) d ^ for any % G {1, . . . , N}, where > depends on the dimen- 
sion d only. 

(hi) Let v be soft-core, assume that d £ {2, 3}, and let (hi, . . . , h^) be any 
minimizer. Assume that v\(o,n) — f or some 7] > 0. In d = 3, furthermore 
assume that 

(1.13) / \v (\y\)\ 1+S dy < oo for some 5 > 0. 

■ABi(O) 

Then every hi is positive everywhere in M. d and continuously differentiable, 
and all first partial derivatives are a- Holder continuous for any a < 1. 

(iv) Let v be hard-core, assume that tie {2,3}, and let (hi, . . . , /tjv) be 
any minimizer. Then every hi is continuously differentiable in the interior 
of its support, and all first partial derivatives are a-Holder continuous for 
any a < 1 . 

For the proof see Sections 4.1-4.4. 

Remark 1.4. (i) Unlike for the ground states of Hn in (1.6), there 
is no convexity argument available for the formula in (1.9). This is due 
to the fact that a convex combination of tensor-products of functions is 
not tensor-product in general, and hence, the domain of the infimum in 
(1.9) is not a convex subset of H 1 (M. dN ). However, for /12, . . . , hjy fixed, the 
minimization over hi enjoys the analogous convexity properties on i/ 1 (!R' i ) 
as the minimization in (1.6). 
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(ii) If v is hard-core, it is easy to see that the distances between the 
supports of hi, . . . ,/itv have to be no smaller than a [see (1.2)] in order to 
make the value of (hi <g> • • • <g> /i^r, Wtv^i (8) • • • ® hjy) finite. The potential 
J2j^i Vh"j is equal to oo in the a-neighborhood of the union of the supports 
of hj with j / i, and hi is equal to zero there (we regard • oo as 0). In 
particular, minimizers of (1.9) are not of the form (h, . . . , h). This shows that 
the inequality in (1.10) is even strict. In the soft-core case, this statement 
is not obvious at all. A partial result on this question in d = 3 will be a 
by-product of Section 1.5 below. 

1.3. The Brownian models. We introduce two different models of inter- 
acting Brownian motions. These models are given in terms of transformed 
measures for paths of length (3 in terms of certain Hamiltonians. Let a family 
of N independent Brownian motions, (B^ 1 )t>o, ■ ■ ■ , (b\ )t>o> in ^ d with 
generator —A be given. The Hamiltonians of both models possess a trap 
part and a pair-interaction part. The trap part is for both models the same, 
namely, 

N rP 

(1.14) H N> p = J2 W(B^)ds. 



We assume that the joint starting distribution of the motions, that is, the 
distribution 
subset of fL 



distribution of the vector (B^\ . . . , Bq N ^), is concentrated on a compact 



1.3.1. The canonical ensemble model. The Hamiltonian of our first model 
consists of two parts: the trap part given in (1.14), and a pair-interaction 
part, 

(1.15) G NJ3 = J2 fv(\B^-B^\)ds. 

l<i<j<N J ° 

We look at the distribution of the N Brownian motions under the trans- 
formed path measure 

t& NtP = — — exp(-HN,p ~ G N ^) dP, 

(1.16) 

where Z Nj p = E(exp(-H Nj p - G Ni p)). 

We call Pat/3 the canonical ensemble model, since it is derived, via a Feynman- 
Kac formula, from the trace-class operator of the canonical ensemble, e~^' HN 
(see Remark 1.6 below). It is a model for N Brownian motions in a trap 
W with the presence of a repellent pair interaction. We can conceive the 
iV-tuple of the motions, B t = (S t (1) , . . . , B { t N) ), as one Brownian motion in 
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M. dN . Introduce the normalized occupation measure of the diV-dimensional 
motion, 

1 [f 3 

(1.17) np(dx) = - 5 Ba {dx)ds, 

P Jo 



which is a random element of the set .A/fi(R ) of probability measures on 
M. dN . It measures the time spent by the tuple of N Brownian motions in 
a given region. Note that there is only one time scale involved for all the 
motions, that is, the Brownian particles interact with each other at com- 
mon time units. We can write the Hamiltonians in terms of the occupation 
measure as 

(1.18) H Nl p = P{im,ni3) and G Ny g = /3(t>, up), 

where the functions 2H, d :R — ► R are introduced in (1.3). 

It turns out that the large-/? behavior of the canonical ensemble model is 
described by the ground state of the operator TLn in (1.7). In Section 1.6.2 
below we recall the notion of a principle of large deviations. The rate func- 
tion Jjv appearing in Theorem 1.5 is the well-known Donsker-Varadhan rate 
function on M. dN defined by 



(1.19) Jjv(/i) 




if ./W eH 1 (R dN ) exists, 
V ax 

otherwise. 



Note that the energy functional (hjTijfh) may be rewritten (h, Ti^h) = 
In(i^) + (2B? A 4 ) + ( D i A*) f° r the probability measure fi(dx) = h?{x) dx. 

Theorem 1.5 (Canonical ensemble model at late times). Fix JV6N. 
(i) 

(1-20) lim -L logE(ex P (-F JVl /3 - G N j,)) = - X N, 

where xn is the ground- state energy per particle of the N -particle operator 
Hn given in (1.6). 

(ii) As (3— >oo, the distribution of fj,p on 7Wi(R dJV ) under F N> p satisfies 
a principle of large deviation with speed (3 and rate function In given by 

(1.21) I N (it)=I N (j J ) + (W,ii) + (x>,ii)-NxN for At € Mi (R dAr ). 

(iii) The distribution of fig under Pjv g converges weakly toward the mea- 
sure h*(x) 2 dx, where /i* is the unique minimizer in (1.6). 



For the proof see Section 5. 
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Remark 1.6. It is well known [18] that the bottom of the spectrum of 
TCn is related to the large-/? behavior of the trace of e~^ N , more precisely, 

(1.22) X N = - Jim -^logTV(e-^). 
Using the Feynman-Kac formula for traces [7, 18], we have 

(1.23) Tr(e"^)= / dxE^ ( e -W->M) ; 

Jn 

where is the expectation with respect to a Brownian bridge of length 
(3 that starts and ends at x S M. dN . Note the close relation to Theorem 1.5. 

1.3.2. The Hartree model. Our second Brownian model is defined in 
terms of another Hamiltonian. We keep the trap Hamiltonian H^^ as in 
(1.14), but the interaction Hamiltonian is now 

(1.24) K N ,p= J2 \ff P v(\Bf)-B[ j) \)dsdt. 

Note that the ith Brownian motion interacts with the mean of the whole 
path of the j'th motion, taken over all times before (3. Hence, the interaction 
is not a particle interaction, but a path interaction. The interaction (1.24) 
is related to Polaron-type models [5, 13], where instead of several paths a 
single path is considered. We consider the corresponding transformed path 
measure, 

<*S?J = -gj exp(-H,,,p - A>,„) <(P, 

where Z y N £ = E(exp(-i/ 7V]/3 - K N ^)). 

In Theorem 1.7 below it turns out that the large-/? behavior of Z^p is 
intimately related to the Hartree formula in (1.9). Therefore, we call this 
model the Hartree model. At the end of this section we comment on its 
physical relevance. 

We introduce the normalized occupation measure of the ith motion, 

.<«),. • 1 fP 



(1.26) ^{dx) = -J^ 5 B w(dx)dseMi(M d ). 

The tuple of the occupation measures, {/J>g \ ■ ■ ■ ,/A ), plays a particular 
role in this model. We can write the Hamiltonians as 

N 

HN,f} = P52(W,iif)=l3{<nJ,ij!$) and 

(1.27) i=i 

l<i<j<N 
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where we recall (1.3) and the operator V with kernel v o | • |, and //® 

/ui <8> ■ • ■ <S> ftp is the product measure. Now we formulate the main result 
on the large-time behavior of the Hartree model. 



Theorem 1.7 (Hartree model at late times). Assume that W and v are 
continuous in {W < oo}, respectively, in {v < oo}. Furthermore, assume 
in the soft-core case that there exists an e > and a decreasing function 
v:(0,e) with v <v on (0,e), which satisfies f B j Q \ G(0,y)v(\y\) dy < oo, 

where G denotes the Green's function of the free Brownian motion on M. d . 
Fix N G N. 

(i) For as defined in (1.9), 
(1-28) Jim ^l og E(e X p(-F iV>/3 - K NtP )) = - X f ] - 

As (3 — > oo, the distribution of the tuple (fJa , • • • , Ma ) of Brownian 

occupation measures on A4i(M. d ) N under satisfies a large deviation 

principle with speed (3 and rate function 



N 

if ] (m,...,fi N ) = J2h Oh) + (an, m®) + <* , m®) - n Xn , 



(1-29) i_1 

fi u ...,fi N £M 1 (M. d ), 

where I\ is defined in (1.19), and fjP = fj,i ® • • • ® /xjv is the product measure. 

(hi) The distribution of (pSp , . . . , /A , ) under is attracted by the set 
of minimizers in (1.11). 



For the proof see Section 2. 



Remark 1.8. (1) The additional assumption of continuity of W and v 
is necessary only in the proof of the upper bound, where we rely on large- 
deviation arguments and need continuity of the map fi ([i,W A M + (v o 
| • |) A M) in the weak topology on the set of probability measures. In the 
proof of the lower bound, we use an eigenvalue expansion, which needs only 
local integrability. 

(2) The additional assumption that there is a function v with v < v and 
/bi(o) ^(0> v )v{\y\) dy < oo in the soft-core case is necessary only in our proof 
of the lower bound. This means that we can handle the case where v(r) < 
0(r~ e ) as r | with e < 2, but not the case v(r) > Cr~ e with e E [2,d], as 
the sole requirement f Bl ^ v(\y\) dy < oo would include. 
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Let us draw a corollary about the mean of the N Brownian occupation 
measures, 



(1.30) 



i N 

1 (i) 



i=l 



Corollary 1.9. Fix iV e N. As (3 — ► oo, the distribution of the mean 
/Ijv,/3 under Pj^l satisfies a large deviation principle in the weak topology on 
Mi(M. d ) with speed j3N and rate function 

N 



(1.31) 



+ inf 

)j,l,...,fi N eMi{R d ): li=IJ. 



In particular, /ijv,/3 converges under Pj^g, as (3 — > oo, weakly toward the 
minimizer(s) o//j\f' mean ^ on Mi(M. d ). 



N 



N 



i=l 



N 



2N 1 



i=l 



(®) 



Remark 1.10. If one changed the Hartree model by adding all the terms 
for i = j in the pair-interaction term, then all statements would remain valid 
after obvious changes in the notation. This model has an additional self- 
interaction of each path. For this model, stronger statements are possible. 
Assume that v is such that the map fiH (/x,V//) is convex on Ai\(M. d ). 
Then the rate function in Corollary 1.9 (without the term V^i)) 
may be identified as I(fi) + (W, fi) + (/i, Vfi) (minus the normalization) and 
therefore turns out to be strictly convex. 



Let us now comment on the physical relevance of the Hartree model. 
Recall that the Hartree formula in (1.9) was introduced as an ansatz for 
studying the ground states oi TLn- Theorem 1.7 shows that the Hartree 
model is the correct positive-temperature model for this ansatz. In this way, 
its relation to the product ground states of Wat is analogous to the relation 
of the canonical ensemble model to its ground state. 

The study of the large- N behavior of the canonical ensemble model at 
positive temperature is an important and difficult open problem. We con- 
ceive the Hartree model as an ansatz for approaching this task. Indeed, we 
study this question for the Hartree model in [1]. For this question at zero 
temperature, see Section 1.5. 



TRAPPED INTERACTING BROWNIAN MOTIONS 



13 



1.4. Dirac interaction. In this section we discuss a further interesting 
choice of the interaction potential. We restrict now to dimensions d E {2,3}. 
We do not choose the interaction as a function, but as a measure, still keeping 
the singularity at zero. More precisely, we replace v by the Dirac measure 
at zero, 5q. 

Let us first remark that the canonical ensemble model does not feel this 
interaction, since P(3 s:Bi^ = Bg) = and hence, J2i<i<j<N fo $o(\Bs — 

Bs \)ds = almost surely. Hence, this model with v replaced by 5o is the 
same as if v would be replaced by 0. 

However, in d = {2,3}, the Hartree model is nontrivial and highly inter- 
esting. In order to see this, recall the intersection local time of the ith. and 
the jth motion for i < j, which is formally defined by 

(1.32) af\x) = ±^ds J*dt8*{\B®-BP\). 

It is known [17] that there is a continuous stochastic process (<Xa { x )) x eS. d 
which justifies this formal definition, that is, for any continuous bounded 
function / — > M, we have the formula 

/ f(x)af j \x)dx = ± fds f dtf(\Bf-B^\). 

JB. d H p JO JO 

(i i) 

Furthermore, a a is even continuous in 0, which makes it possible to de- 
fine the normalized amount of intersection of B^ and B^ as the random 
variable al' J '(0). Hence, we can replace the Hamiltonian Kjy^ in (1.25) by 

(1.33) K N ^ = X(3 J2 a f j \°) for 

some A £ (0, oo). 

l<i<j<N 

The relevant variational problem should be 

f N 

/ -i a\ x inf < > llV/iilU 

l i - d4 J h 1 ,...,h N £H 1 (R d ): ||hi|| 2 =lVi {~[ 

1 N 1 

i=l l<i<j<N ) 

Using the means of the proof of Lemma 1.3 in Sections 4.1-4.3 below, one 
can also show the following. 

Lemma 1.11 [Minimizers in (1.34)]. Fix A 6 (0, oo). 
(i) The infimum in (1.34) is attained. 
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(ii) Any minimizer (h\, . . . ,hpf) of (1-34) satisfies the following system 
of Euler-Lagrange equations: 

Ahi = -\ihi + Whi + Xhi^hj, i = l,...,N, 

where Xi = \\Vhig + {W,h}) + X E^i(^f, h)) . 

(iii) Let (h±, . . . , /tj\r) be a minimizer of (1.34). Then every hi is positive 
almost everywhere. Furthermore, in d = 2, every hi is continuously differ- 
entiable, and the first derivatives are a-Holder continuous for any a < 1. In 
d = 3, every h\ is locally integrable for any p 6 (0, 3). 

It is natural to conjecture that, for any N € N and any A € (0, oo), 

(1.35) lim J^logE( e -^-^) = -x!J o) (A), 

(B— >oo I\ p 

and that the tuple of normalized occupation measures of the motions stands 
in the analogous relation to the minimizers as in Theorem 1.7. However, we 
do not know how to prove this conjecture, since the treatment of Brownian 
intersection local times is technically notoriously difficult. Instead, we offer 
an analogous result for simple random walks in place of Brownian motions, 
which we can handle rigorously. 

For this purpose, let )te[o,oo) be independent continuous-time simple 
random walks on T, d for i = 1, . . . , N with generator —A given by Af(z) = 
J2y~z(f( v ) ~ f( z ))- For simplicity, we pick some joint initial distribution 
with compact support. 

The normalized intersection local time of the ith and the jth walk is given 

as 

where £$(z) = Jq 1{S^ = z}dt denotes the local times of the ith walk up 
to time (3 > 0. In order to have a perfect-scaling property, we restrict to the 
trap potential W(x) = \x\ p for some p > d — 2. We consider the following 
model: 

1 / 1 N rP 
dipW =-(^exp --W W(S®)dt 

(1.36) Z ^ V P *= l ° 

_ X(3 (d +P )/(2 +P ) af j) \d¥. 

l<i<j<N ) 

As usual, Zjya > denotes the constant that makes \ a probability mea- 
sure. Note the factor of 1/(3 in front of the trap term and the factor of 
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/3(d+p)/( 2 +p) in front of the interaction term. The first ensures that the prop- 
erly rescaled random walks approach some Brownian motions, which makes 
the model asymptotically equal to the above Brownian model. In order to 
formulate our result on the large-/? asymptotic of P$g, we nee d to introduce 

(i) 

the following normalized and rescaled version of 1^ . 

(1.37) Lf{x) = ^lf{\x^\), x e R d , where & = p 1 /^. 

Note that Lf is a (random) probability density on W 1 . Let /Jg\dx) = 
L$ (x) dx be the corresponding measure. 

Theorem 1.12 (Discrete Dirac-interaction model at late times). Fix 
d>2 and p > d — 2. Furthermore, let N G N and A £ (0, oo) . Then 



(i) 

(1.38) lim /T p/{2+p) log 4^ = -Nx { n ] W- 



As (3^oo, the distribution of the tuple (/A , • • • , Ma ) of normal- 
ized rescaled local times under on Mi(M. d ) N satisfies a large deviation 
principle with speed (3 P ^ ( 2+p ) and rate function 



N 



I^^ u ..., m ) = J2h(^) + (W,^) 
i=i 

(1-39) +A £ (f^)-A#(A), 

^,...,^ N eMi(R d ), 

where I\ is defined in (1.19), )U® = fj,± <8) • ■ • <8) //at «s i/ie product measure, 
and we define . . . , /ijv) = +oo z/ any o/ i/ie measures . .. ,/ijy /a^s 

to /iai>e a Lebesgue density. 

(iii) T/ie distribution of (/A , • • • , under « is attracted by the set 

of minimizers in (1.34). 

For the proof see Section 3. 

Note that Theorem 1.12 is true in any dimension d>2, while the Brow- 
nian version is well defined only in d £ {2, 3}. 

The choice of the factor 1/(3 in front of the trap term in (1.36) is high- 
handed. Because of this term, one has to assume that p > d — 2 in order that 
an appropriate large deviation principle be applicable (see Lemma 1.16). If 
1//3 would be replaced by (3~ 2 ^ d , then the assumption p > would suffice. 
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1.5. Large-N behavior of the product ground states. In this section we 
study our main variational formulas, xn and xjy\ an< ^ their minimizers in 
the limit for diverging number N of particles. In particular, we point out 
some significant differences between \N and its product state version yffi 
in the soft-core and the hard-core case, respectively. 

First we report on recent results by Lieb, Seiringer and Yngvason on the 
large- N behavior of xn- Let the pair functional v be as in (1.2) and assume 
additionally that v > and v(0) > 0. 

We shall replace v by the rescaling vn(-) = £jf 2 v(- £jv )> f° r some appropri- 
ate £jv tending to zero sufficiently fast. Hence, the reach of the repulsion is 
of order £tv, and its strength of order £^ 2 . Furthermore, the scattering length 
of v, a(v), is rescaled such that a(u7v) = cz{v)(3n (see Section 1.6.1 below for 
the definition of the scattering length and some of its properties). If (3m I 
sufficiently fast, this rescaling makes the system dilute, in the sense that 
a(v N ) < AT-Vf This means that the interparticle distance is much bigger 
than the range of the interaction potential strength. More precisely, the de- 
cay of Pn will be chosen in such a way that the pair-interaction has the same 
order as the kinetic term. In d = 3, this choice of rescaling is motivated by 
famous experiments for the derivation of Bose-Einstein condensation of a 
large, but finite, dilute trapped system of N real particles ( 87 Rb [4], 7 Li [6], 
23 Na [10], but also more recently 85 Rb, 41 K, 133 Cs, hydrogen, metastable 
triplet 4 He, 174 Yb, 85 Rb2, and 6 Li2). Here the scattering length is of the 
order 10 -3 , whereas N varies from 10 3 to 10 7 . 

The mathematical description of the large- iV behavior of xn in this scal- 
ing, and hence, the theoretical foundation of the above mentioned physical 
experiments, has been successfully accomplished in a recent series of pa- 
pers [22, 24, 25, 26]. It turned out that the well-known Gross-Pitaevskii 
formula adequately describes the limit of the ground states and its energy. 
This variational formula was first introduced in [19] and [20] and indepen- 
dently in [27] for the study of superfluid Helium. After its importance for 
the description of Bose-Einstein condensation of dilute gases in magnetic 
traps was realized, the interest in this formula considerably increased; see 
[9] for a summary and the monograph [28] for a comprehensive account on 
Bose-Einstein condensation. 

The Gross-Pitaevskii formula has a parameter a > and is defined as 
follows: 

(1.40) xi GP) = inf [||V^|| 2 + (W,^)+47ra||^|| 4 ]. 

< ? iGJ/ 1 (M d ) : ||0|| 2 = 1 

It is known [22] that Xa* possesses a unique minimizer 4>^ P \ which is 
positive and continuously differentiable with Holder continuous derivatives 
of order one. 
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Since v(Q) > 0, its scattering length a(v) is positive (see Section 1.6.1 
below). The condition f^ l v(r)r d ~ 1 dr < oo implies that a(v) < oo. Fur- 
thermore, note that the rescaled potential £~ 2 v(- 1;" 1 ) nas scattering length 
£ot(v) for any £ > 0. 

Theorem 1.13 (Large-iV asymptotic of \N in d G {2,3}, [22, 24, 26]). 
Assume that d G {2, 3}, i/iaf t> > with v(0) > 0, and J^_ 1 i)(r)r d ~ 1 dr < oo. 
Replace v by ujv(-) = £,n 2 v(- ^ 1 ) wzi/i £tv = 1/iV in d = 3 and £jy = a(w)~ 2 x 
g -A r /«(j;) jy|| ( ^ ) (G ; P)||-4 ^ n ^ — 2. £ e £ /ijy G i7 1 (IR ciAr ) 6e £/je unique minimizer 

on the right-hand side of (1.6), and define (j) N G i? 1 (]R d ) as i/ie normalized 
first marginal of h 2 N , that is, 

4>%(x)= h 2 N (x,X2, ■ ■ ■ ,XN)dx 2 - ■ -dx N , x£R d . 

J R d(N-l) 

Then we have 

lim XN = Xa(v) and $ 2 N^(4>aMf in weak L l {W 1 )- sense. 

In particular, the proofs show that the ground state, Hn, approaches 
the ground state {4>^^)® N if N gets large. In order to obtain the Gross- 
Pitaevskii formula as the limit of xn also in d = 2, the rescaling of v in 
Theorem 1.13 has to be chosen in such a way that the repulsion strength 
is the inverse square of the repulsion reach and such that this reach decays 
exponentially, which is rather unphysical. 

In the present paper, we prove the analogue of Theorem 1.13 for the 

Hartree model in the soft-core case. It turns out that Xn m (l-H) a l so 
converges toward the Gross-Pitaevskii formula. However, in d = 2, it turns 
out that the potential v has to be rescaled differently. Furthermore, in d G 
{2,3}, the scattering length a(v) is replaced by the number 

(1.41) 5(«):=-!- / v(\y\)dy. 

Theorem 1.14 (Large- N asymptotic of Xn ' soft-core case). Let dG 
{2,3}. Assume that v is a soft-core pair potential with v>0 and v(0) > 
and a(v) < oo. In dimension d = 3, additionally assume that (1.13) holds. 

Replace v by ujv(-) = N d ~ 1 v(- N) and let (h± , . . . , fiff ) be any minimizer 

on the right-hand side of (1.11). Define 4> 2 N = i J2iLi(h\ N ^ ) 2 ■ Then we have 

lim xn =X^\ an d 4>n ~ * ) 2 > 

where the convergence of (f> N is in the weak L l (M. d ) -sense and weakly for the 
probability measures (j) 2 l (x)dx toward the measure (c/>~^) 2 (x) dx. 
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For the proof see Section 4.5. 

Note that, in d = 3, the interaction potential is rescaled in the same way 
in Theorems 1.13 and 1.14. However, the two relevant parameters depend 
on different properties of the potential (the scattering length, resp. the inte- 
gral) and have different values, since a(v) < a(v) (see [26] and Section 1.6.1 
below). In particular, for N large enough, the ground state of xn is not a 
product state. This implies the strictness of the inequality in (1.10), for v 
replaced by vn(-) = N 2 v{- N). The phenomenon that (unrestricted) ground 
states are linked with the scattering length has been theoretically predicted 
for more general iV-body problems [14], Chapter 14, [29]. Indeed, Landau 
combined a diagrammatic method (a Born approximation of the scattering 
length) with Bogoliubov's approximations to almost reconstruct the scat- 
tering length from the L 1 -norm of v o | • | in the (nondilute) ground state. 
However, the relation between the L 1 -norm and the product ground states 
was not rigorously known before. 

In d = 2, a more substantial difference between the large-iV behaviors 
of xn an d Xn is apparent. Not only the asymptotic relation between the 
reach and the strength of the repulsion is different, but also the order of this 
rescaling in dependence on N. We can offer no intuitive explanation for this. 

Interestingly, in the hard-core case, Xm shows a rather different large- 
N behavior, which we want to roughly indicate in a special case. Assume 
that W and v are purely hard-core potentials, for definiteness, we take W = 
°°1bi(o) c an d v = ool[o )0 ]. We replace v by uat(-) = v{-/£,n) for some £jv I 

(a pre- factor plays no role) . Then x%^ is equal to times the minimum over 
the sum of the principal Dirichlet eigenvalues of — A in N subsets of the unit 
ball having distance > cl(3n to each other, where the minimum is taken over 
the N sets. It is clear that the volumes of these ./V sets should be of order 
jj, independently of the choice of Then their eigenvalues are at least of 

order N 2 / d . Hence, one arrives at the statement liminfAr_>oo N~ 2 l d x^ > 0, 

that is, Xtv^ tends to oo at least like N 2 / d . 

1.6. Preliminaries. 



1.6.1. The scattering length. Let us briefly introduce the scattering length 
of the pair potential, v, and its most important properties. For a detailed 
overview, see [26]. First we turn to d > 3. Let u: [0, oo) — ► [0, oo) be a solution 
of the scattering equation, 

„," _ l, 



(1.42) 



juv on (0, oo), u(0) = 0. 



Then the scattering length a(v) £ [0, oo] of v is defined as 

u(r) 



(1.43) 



a(v) 



lim 

r— >oo 



u'(r) 
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If v(0) > 0, then a(v) > 0, and if J^_ 1 v(r)r d ~ l dr < oo, then a(v) < oo. In 
the pure hard-core case, that is, v = ool[ . a ), we have a(v) = a. It is easily 
seen from the definition that the scattering length of the rescaled potential 
£ _2 w(-£ _1 ) is equal to £a(v), for any £ > 0. 

There is some ambiguity of the choice of u in (1.42); positive multiples of 
u are also solutions, but the factor drops out in (1.43). We like to normalize u 
by requiring that liniR_ >00 ?/(.R) = 1. It is easily seen that (where oJd denotes 
the area of the unit sphere in 



u (W)i — \lZo^ x = UJ d I v{r)u{r)r dr 







oo 



= 2u>d / u (r)rdr 
(1.44) d Jo U 

= 2uj(i lim iu(r)r\Q— / u'(r)dr 
R^co \ Jo 

= 2uj d lim (u(R)R-u(R)) = 2u> d a{v). 

R— too 

As a consequence, in dimension d = 3, we have a(v) < a(v). Indeed, u is 
a nonnegative convex function whose slope is always below one because of 
limR^oo u'{R) = 1. By u(0) = 0, we have that u(r) < r = r d ~ 2 for any r > 0. 
With the help of (1.44), we therefore get 87ra(f ) = 2ujda{v ) < J Rd v(\x\)dx = 
8ira(v). 

In d = 2, the definition of the scattering length is slightly different. We 
treat first the case that supp(u) C [0,i2*] for some iZ* > and consider, for 
some R> R*, the solution u: [0,R] — > [0,oo) of the scattering equation 

u" = \uv on[0,R], u(R) = l,u(0)=0. 

Then u(r) = log / log for R* <r < R for some a(v) > 0, which is by 

definition the scattering length of v in the case that supp(w) C [0,iZ*]. Note 
that a(v) does not depend on R. Hence, 

logr -u(r) log R 

\oga(v) = — , R*<r<R. 

1 — u{R) 

For general v (i.e., not necessarily having finite support), v is approximated 
by compactly supported potentials, and the scattering length of v is put 
equal to the limit of the scattering lengths of the approximations. 

1.6.2. Large deviations principles. For the convenience of our reader, we 
repeat the notion of a large deviation principle. A family (X^) / g > g of random 
variables Xp, taking values in a topological vector space X, satisfies the large 
deviation upper bound with speed ap, where ap — > oo for (3 — > oo, and rate 
function I : X — > [0, oo] if, for any closed subset F of X, 

lim sup — log¥(Xp E F) <— inf I(x), 

/3->oo ap xGF 



20 



S. ADAMS, J.-B. BRU AND W. KONIG 



and it satisfies the large deviation lower bound if, for any open subset G of 
X, 

liminf — logP(Xg GG)>- inf I(x). 

/3— >oo dp x&G 

If both, upper and lower bound, are satisfied and, in addition, the level sets 
{/ < c} are compact for any c£K, then one says that (Xp)p satisfies a 
large deviation principle. This notion easily extends to the situation where 
the distribution of Xp is not normalized, but a sub-probability distribution 
only. 

In the proofs of Theorem 1.7 we shall rely on the following principles for 
the normalized Brownian occupation measures, that is, for certain .Mi (Re- 
valued random variables. For any measurable subset A of M. d , we conceive 
M\ (A) as a closed convex subset of the space A4(A) of all finite signed 
Borel measures on A, which is a topological Hausdorff vector space whose 
topology is induced by the set C\,(A) of all continuous bounded functions 
A -> R. Here C h (A) is the topological dual of M(A). The set M±(R d ) inherits 
this topology from Ai(A). When we speak of a large deviation principle of 
Mi(^4)-valued random variables, then we mean a principle on M(A) with a 
rate function that is tacitly extended from M\{A) to M(A) with the value 
+00. 

One of the principles we are going to present is for the Brownian motion 
in a given bounded set, and the other one for a periodized version of the 
motion, which we introduce now. Let R > and A R = [-R, R] d C R d . For 
any probability measure \x on U. d , we denote by hr € Mi(Ar) the periodized 
version of //, that is, 

(1.45) Hr{A) = //I (J (A + 2kR) J , A C A R measurable. 

Note that the shifted cubes Ar + 2kR with k G 7L d are disjoint up to their 
boundaries, and that their union covers M. d . Recall the Donsker-Varadhan 
rate function from (1.19). The periodized version of the rate function I\ on 

M\{Ar) is denoted , that is, 

(1.46) /rV) = mf{/iH : /i = ur is the periodized version of v}. 

Lemma 1.15 (Large deviations principle for occupation measures [16]). 
Fix d G N. Let {Bt)t>o be the Brownian motion on M. d with generator —A, 
and let fip(dx) = 4 Jq ' 5e s (dx) ds be the normalized occupation measure up 
to time (3 > 0. Fix R>0. 

(i) The family (/^ / 3 ! _r)/j>o of An-periodizations of [x^ satisfies the large 
deviations upper bound with speed (3 on A^i(A^) with rate function 1^ . 
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(ii) For any open bounded set A C M. d , the family (iip)/3>o satisfies, un- 
der the sub-probability measures P( • n{supp(/ig) C A}), the large deviations 
lower bound with speed (3 on Mi (A) ; the corresponding rate function is the 
restriction of I to the set of probability measures whose support lies in A. 

In the proof of Theorem 1.12 in Section 3 we shall rely on related princi- 
ples for the normalized and rescaled local times defined in (1.37), more 

precisely, for the corresponding measures fj?p\dx) = L^\x) dx. 

Lemma 1.16 (Large deviations principle for rescaled local times of ran- 
dom walks). Fix d € N. Let (Sj)te[o,oo) be a simple random walk on 7j d 
with generator —A (the discrete version of the Laplace operator) and local 
times lp(z) = Jq t{S t = z] dt, and let 1 -C £p <C P l ^ d as (3 — > oo be some scale 
function. In d = 1 assume that £p <C \J~P, in d = 2 assume that p/ 'log p. 
Define the normalized and rescaled occupation measure by 

W(dx) = -£i(\x£p\)dx. 

(i) The family (/x ( g ) i?) l g>o of An-periodizations of [ip satisfies the large 
deviations upper bound with speed Pip 2 on Mi(Ar) with rate function l[ R \ 
(ii) For any open bounded set A C M. d , the family (ug)/j>o satisfies, un- 
der the sub-probability measures P(- D {supp(/ig) C A}), the large deviations 
lower bound with speed P£p 2 on Mi(A); the corresponding rate function is 
the restriction of I to the set of probability measures whose support lies in 
A. 

Proof. See [15], Lemma 3.2, for the case of a discrete-time random 
walk. The proof for the continuous-time setting is very similar and is there- 
fore omitted. □ 

2. Large deviations for the Hartree model: Proof of Theorem 1.7. In 

this section we prove Theorem 1.7. We shall proceed according to the well- 
known Gartner-Ellis theorem. Therefore, we have to establish the existence 
of the logarithmic moment generating function of (/A 1 , . . . ), that is, 
the existence of 

(2.1) A<f >(*) = foa i logE^exp^GuW . . .J^))], 

for any element <5 of the dual of the vector space M(M. d ) N . Note that ev- 
ery linear continuous functional on Mi(M. d ) N is of the form (/xi, . . . ,/x/v) l— ► 
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12i=i(fii Pi) with /i,...,/jv £Cb(M d ), the set of bounded continuous func- 
tions on R rf . 

The main step in the proof of Theorem 1.7 is the following. 

Proposition 2.1 (Asymptotic for the cumulant generating function). 
For any h,...,f N er- 



(2.2) 

where 

(2-3) X f\f) 



1 



lim -logE[e _H 

/3— »oo p 



NxfHf), 



— inf 

W fn,...,^ N eMi 



N 



^I 1 (/x i ) + (2H + D-/,^ 



i=l 



and we wrote //® = /ii <8> • • • <8> ^tv / = /i © • • • © /jv- 



The next three subsections are devoted to the proof of the upper and 
lower bound in (2.2), respectively. In Section 2.4 we finish the proof of The- 
orem 1.7. 

An outline of our proof is the following. Recall (1.26) and (1-27) to see 
that 



E\e 



-H 



N,0 



-K 



(2.4) 



exp< 



N 

1=1 



ft 



E 

l<i<j<N 



i=l 



(*)\ 



:E[exp{-/3(2U+o-/, M |)}] 



We intend to apply the large deviation principle in Lemma 1.15 and Varad- 
han's lemma, which would immediately yield the result in (2.2). However, 
there are some technical obstacles to be removed. Because W explodes at in- 
finity and v has a singularity at the origin, both functionals ii/jv/3 an d Kn,/3 
are not continuous and not bounded in the weak topology. Furthermore, we 

cannot apply immediately the large deviation principle of Lemma 1.15, be- 

(i) 

cause the occupation measures /A are not restricted to any bounded domain 

in M. d . As it concerns the proof of the upper bound, these technical obstacles 
will be removed in Section 2.1 via a well-known cutting and periodization 
procedure using the large deviations principle in Lemma 1.15 (i) . An analo- 
gous technique works for the proof of the lower bound in the hard-core case, 
using Lemma 1.15(h). However, for proving the lower bound in the soft-core 
case, we did not succeed in making the principle in Lemma 1.15(h) applica- 
ble. The main reason is the singularity of v at zero, which seems to destroy 
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all necessary semicontinuity properties. Instead, we employ an eigenvalue 
expansion technique for the iV iterated expectations w.r.t. the N motions. 
Here the additional integrability property of v is necessary. 

2.1. Proof of the upper bound in (2.2). We consider the large closed box 
Ar = [-R, R] d . We divide the probability space into the part on which each 
motion spends more than (1 — rj)[3 time units in Ar up to time (3 (the main 
part) and the remaining part, where this is not satisfied (this part will turn 

out to be negligible). On the first part, we shall replace each //g by its 

A^-periodized version fJ$ R and control the error. Then -ff/v,/? turns out to 
be a continuous and bounded functional of the periodized versions. Also, 
we replace the functions W and v by their cut-off versions Wm = W A M 
and vm = v A M, respectively, where M > is large. This will enable us to 
apply Varadhan's lemma. Finally, we let the auxiliary parameters R, r\ and 
M tend to infinity respectively to 0. 

We turn to the details. Let auxiliary parameters rj > be small and 
R, M > be large. The expectation on the left-hand side of (2.2) is split 
into two parts and get the estimation 



^ e -H Ni p-K N>p +/3{f,fip v 

(2.5) < E ( e -^-^W^?) Y[ l{fif (A R ) > 1 - r?} 



N 



i=l 



+ NE(e- H »>n{$\A R ) < 1 - ?? }) e -/3[(^/2)mf,-c / ] ) 



where Cf = Ya=x WfiWoo- The second term is easily estimated, using that 
W > 0. Indeed, we have, on the event {^'(Ar) < 1 — rj}, 

H N ,f3 >p f W{x l )^\dx 1 ) > pqipfW, 
Ja c r a r 



and therefore, the second term on the right-hand side of (2.5) is not bigger 
than Ne -to™# a W-W*/*)***-C,)_ 

In the first term, we first estimate K^,p = /3(t),//®) > P(0m,^}, where 
vm{u) = v(y) A M is the cut-off pair potential, and Dm is defined as t> with 
v replaced by vm- Analogously, we estimate = (3(W, ^} > /3(Wm, ^) 
with analogous notation. This leads to 

^ e -H Ntl 3-K Ntl3 +f3{f,Li®) j 

< Ne -P(r,ird A c R W+(N 2 /2)infv-C f ) 

(2-6) " N 

+ E L-W»,$)-M*M,t$)+W$) ftt{$(A R ) >l-v}Y 
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Now we replace by its periodized version /if R = l^p R ® • ■ • ® Mg In 
order to estimate the error, we point out that, for any i £ {1, . . . ,N}, 

/ W M {Xi)Hg\dXi) - W M {Xi)^a R {dXi) 

Jm d Ja r 

= J2 I (W M (x i + 2Rk)fif(d(x i + 2Rk)) 

{ > -W M {xi)lif{d(xi + mk))) 

= E / nf(d(xi + 2Rk))(WM(xi + 2Rk)-W M {xi)) 
kez d \{o} Ar 

> -T)M. 

Analogously, we derive the error estimate (f,^®) < The 
replacement error for the second term is estimated in a similar way: for any 
pair i,j £ {1, . . . , N} with i / j, we have 

v M {\xi - Xj\)fJp (dxi)[ij\dxj) 

v M (\xi - Xj\)^ R (dxi)^ R (dxj) 

[v M {\xi-x j +2R(k-l)\)-v M {\xi-x j \)} 
( 2 - 8 ) k,iez d -.k^i jARjAR 

x nf{d{xi + 2Rk))^\d(xj + 2RI)) 
>(Mv-M) J2 vf(^R + 2Rk)Jf{k R + 2Rl) 

k,leZ d : kfX 

> -2rj(M -infu), 

since — oo < infu < vm < M and fi^\A R ) < rj. Summarizing, we obtain 
from (2.6) that 

^ e - H N,f3-K Ntl 3+f3{f,n®)^ 

(2.9) < AT e ~ /3(??infA B W+N 2 /2infv-C f ) 

+ e )3r ] lNM+N 2 (A^-miv)+C f }- E ^ e f3{~W M ~x, M +f,fi^ R )^ 

We now argue that we have 

(2.10) limsupilogE(e /3< - 2I,M - l ' M+/ '^« ) ) < -N X f\R,MJ), 

/3-»oo P 
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where 

M 



(2.11) N 



i r N 

- inf j£ ±||V*M*,2 + ((23 + Dm - /)| A g, (^ ) 2 } : 



^i,...,/uvGC 0O (Afl),||^||fl, 2 = lVi = l,...,iV 



where || • ||fl2 is the norm on L 2 (A#), Vr is the gradient on the torus A# 
(i.e., having periodic boundary condition), and /t® = h\ <g> ■ • ■ (8) ftjv- To prove 
(2.10), we apply the upper bound part of Varadhan's integral lemma; see [11], 
Lemma 4.3.6. Indeed, according to Lemma 1 . 15 (i) , every family (ftp r)/3>o 
satisfies the large deviations upper bound as (3— ► oo, and the map 

Om,.. .,/jln) >->• (-2Um - Dm + (8) • • • <8>Mjv) 

is upper semicontinuous (even continuous) on the set of vectors of probability 
measures on A# in the weak topology. Also, using [11], Example 4.2.7, (2.10) 
is a direct consequence, noting that 



x^(R,MJ) 

(2.12) 



1 f N 

k 2=1 



jf r) (m) + (2H M + Dm - /, in ® ■ • ■ ® A*JV> : 



From (2.9) and (2.10) we obtain that, for any r/,R,M > 0, 

the l.h.s. of (2.2) 

f iV 2 

/ 010 x < — min< 7?inf W H infu — C7, 

(2-13) - \'a^ 2 7 

- V [NM + iV 2 (M - inf v) + C/] + N X f ] (R, M, f) 

On the right-hand side of (2.13), we let first R — > oo, then 77 J, and finally 
M — > 00. Recall that lim^oo inf = 00, according to our assumption in 
(1.1). It is easily seen that the proof of the upper bound in (2.2) is finished 
as soon as we have shown that, for some C > which does not depend on 

M, 

liminfxf \R,M,f) > _!±+ x j®>(M,/) and 
(2.14) M 

hminf X ^ ) (M,/)> X ^ ) (/), 

M— »oo 
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where Xn \M,f) is defined as Xn (/) i n (2-3) with and v replaced by 
Wm and vm, respectively. 

Certainly, we may assume that lim sup^^^ lim sup^^^ Xn ^ f) < 

oo. Introduce 

(2.15) w-X&W)* ' 



inf A ^_ 1 W, 



M 



and note that lim sup^.^^ £r,m < jj for some C > 0, not depending on M. 
To prove the first assertion in (2.14), fix R > and let hi,. . . , £ C 2 (A#) 
be a vector of approximately minimizing functions for the right-hand side 
of (2.11). We shall construct L 2 -normalized functions h\,...,hN G ^ 1 (M a! ) 
such that, for some C > 0, not depending on M nor on R, 



N 



EslWI^ + <(2Ua/ + Dm - /)| A g, (^) 2 ) 
(2.16) i=1 ^ 

> -C^,m + E III v Mi + <2» + °a/ - /, (^ 



i=l 



Passing to the infimum over all vectors of L 2 -normalized functions hi, ■ ■ ■ , /jjv G 
if 1 (]R d ), and letting R — > oo, we then arrive at the first assertion in (2.14). 

In order to show (2.16), pick some (j) € C°°(M. d , [0, 1]) with supp(0) C Ar 
and 0| a_r_i = ^Afl-i- Consider the functions hi = hi<f>/\\hi<f>\\i^2 } trivially ex- 
tended to R d . Hence, hi £ H 1 ^) for any i = 1, . . . , N. Note that 



hj(x) 2 dx < 



J Rd W M (x)h 2 (x)dx 



< — — — Er m- 

mf A c Vkm 



In particular, ||/ij0||2 > 1 — £r,m- 

We now show (2.16). Certainly, we may assume that there is a C > 0, not 
depending on R nor on M, such that |V0| 2 < C and ||V/ii||^ 2 < C fo r an y 
i E {1, . . . , N} and any large R > and any sufficiently large M > 0. Then 
we estimate 
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I^^«IIr,2 — II v^ll! 



IMIIllI^ II 2 



+ 1H^2 f K 1 " tfWH 2 - h 2 \V(f>\ 2 - 2ct>hiVhi • V<j>] 

ll /l i<Pll2 JA-r\Ar-i 



(2.18) - l- e/? , M 

c 



( / h 2 i {x)dx+ f \hi(x)\\Vhi(x)\dx) 

l - £r,m \Ja r \a r - 1 Ja r \a r -i / 



1 — £/?,M 

([ h 2 (x)dx+([ h 2 l (x)dx) \\Vhi\\ 2 ). 

1 — e R,M \JAh\Ar_! \JAr\Ar-! J / 

Now use (2.17) to estimate the right-hand side from below against —£r,mC 
for some C, not depending on M nor on R. 

In order to derive a suitable estimate for the other parts of the functionals, 
note that 

r (i-^ 2 ) (1-HMII5) 



(2.i9) OT a -(o a =E n# n* 



=1 \i=l 



1^112 ||^||2 



l m\2 \\ n m\2 



For the two terms between square brackets in (2.19), use the bounds 

< jj- -772 < iA H \A H _i and 0<— — —2— <- '■ , 

\\hi<t>\\2 Whni l-£fl,Jif 

and recall that lim supjv/_>oo lim sup^^^ Xjv M, /) < oo in order to easily 
arrive at the estimate in (2.16). 

We now prove the second assertion in (2.14). Let (/ii,a/, ■ • ■ , ^7V,Af) be 
minimizers in the definition of xj^(M, /). Along suitable subsequences we 
have, for any i G {1, . . . , iV}, 

^i,M — ► as M — ► oo, 

for some hi G // 1 (R ci ) satisfying ||/ij||2 = 1- The convergence is weak in 
L 2 (R d ), weak for the gradients in L 2 (R. d ), weak in the sense of probabil- 
ity measures and strongly in L 2 (JH d ) on compacts. In particular, 



liminf(i||Vfci,M||! - (f,h 2 M )) > ±||VMi - (f,h 2 ). 

M— »oo 

For any fixed M' > 0, we have 

liminf(/i 2 M ,u M /i 2 M ) > liminf(/i 2 M , v M 'h 2 M ) = (h 2 ,v M 'h 2 ). 

M^oo J M— >oo J 
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Letting M' f oo, the monotone convergence theorem implies that 




In the same way, we see that limmfM-*oo(WM! hf M ) > (W,hf) for any i. If 
we now pass to the infimum over all [hx, . . . ,/itv), we arrive at the second 
assertion of (2.14). This finishes the proof of the upper bound in (2.2). 



2.2. Proof of the lower bound in (2.2), hard-core case. We handle first 
the case when v is a hard-core interaction potential. Recall the parameter 
a = inf{r > 0:v(r) < oo} E [0, oo) from (1.2). Fix a family of open bounded 
sets Ax, . . . , An CM. d whose mutual pairwise distance is bigger than rj, where 
rj > a is close to a. We shall use the lower bound 

^( e - H N,l3-K N ,f3 e l3(f,^)^ 

(2 - 20) > E (><^-/3<^f > e /3(/,Mf > -Q 1{mpp Q i f) c A ^ . 

Note that, on the event n^Li{ su PP(^^) C Ai}, the map 

is bounded since the sets A±, . . . , An are bounded, and since the supports of 
the measures fip , . . . ,fip are bounded away from each other by at least r\ 
[recall from (1.2) that v is bounded on (77,00)]. This map is also bounded and 
continuous on the set of (fix, ■ ■ ■ ,Hn) E A4i(M. d ) N such that supp(/ij) C A; 
for all i. Hence, Varadhan's lemma may be applied and yields, also using 
[11], Example 4.2.7, 

liminf-log[l.h.s. of (2.20)] 

0— >oo p 

( N 

(2.21) > - mf E J (Mi) + <2B + - /, fjP) :n x , . . . , m E Mi(R d ), 

I i=i 

Vi:supp(/ij) C a\. 

Now we substitute hf(x)dx = ^i(dx). According to [8], Proposition 3.29, 
we can restrict to the infimum on hi in C^°(Ai) the set of infinitely often 
differentiable functions M d —> M having compact support in Aj. Also, passing 
to the infimum over all admissible sets Ax, ■ ■ ■ , An, we see that the left-hand 
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side of (2.21) is not smaller than 



1 f N 

-X { N ] (f,v) = "^inf E H V ^ll2 + (9B + * - /, (h 2 



(2.22) 



h l ,...,h N ec™(R d ), 

\\hih = ■■■ = \\h N \\ 2 = l,Vi / j : 
dist(supp(/ij), supp(foj)) > 



where we adapted the notation (h 2 )® = (/ii) 2 <8) • • • <8> (h^) 2 . Now we consider 
the limit as 77 j a. It is clear that 

1 f ^ ^ / \ 

limx^C/, r?) = - inf E II V/n||| + E( W + E ™J " /<> ^ ) : 

(2.23) ^...^jveC^OR"), 

HM2 = • • • = ||/ijv||2 = l,Vi^ j : 

dist(supp(/ij), supp(/ij)) > al. 

Now it is easy to see that the right-hand side of (2.23) is equal to xffi(f)- In- 
deed, let (hi, . . . ,h^) £ i^ 1 (R d ) be a minimizer of the variational formula in 
(2.3). According to Remark 1.4(H) (which is only for fi = 0, but nevertheless 
applies also here), the supports of hi, . . . , hx, which we denote fix, ... , Ojv, 
have distance > a to each other. Furthermore, in the case a = 0, the interi- 
ors of fii, . . . , Oat are disjoint. Approaching hi, . . . , h^ with functions whose 
supports have distances > a to each other, we obtain 



xS?(/) 



I inf (e II vmI + e(^ + E ™? - /i. *?) ■ 

U=l i=l \ / 



/ii £ C c °°(fii) with H/iilla = IVi 1. 

Now one sees that this is equal to the right-hand side of (2.23). This ends 
the proof of the lower bound in (2.2) in the case of a hard-core potential v. 

2.3. Proof of the lower bound in (2.2), soft-core case. Now we turn to 
the proof of the lower bound in (2.2) in the case of a soft-core potential 
v. The proof goes via an iteration of N eigenvalue expansions for the N 
expectations with respect to the N Brownian motions. Let us first handle 
the expectation with respect to the Ath motion, which we denote by . 
We only consider those terms that depend on the Ath motion and work 
almost surely with respect to the first A — 1 motions. 
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First we estimate the expectation under interest in terms of the same 
expectation with respect to some Brownian bridge in the time interval [0, /? + 
1] instead of the Brownian motion on [0,/?]. Let pi(x,y) denote the standard 
transition density of Brownian motion at time 1 from x G M d to y G M d . Note 
that 1 > p\(x,y) for any x,y G M. d . Denote by z^v the initial distribution of 
the iVth motion (i.e., the distribution of Bq N ^) and by the expectation 
with respect to the motion started at x G M. d . We abbreviate q^ = W — fjy + 
T.j<jfVny. Note that (/TV is a random potential which is locally integrable 
in R d , almost surely. Fix some R > 0. Then we can estimate 



> 



v N (dx)EW [e-^^p^p , x)] 
(2.24) > e -\\f N \UHN-D^v f UNidx)E (N )[e-^ +1 ^AVi)s x (B^\)} 

> e -||/iv||oo + (Af-l)inft; 

x / ^(^)4 W) [^ +1)<9 ^ 

The right-hand side of (2.24) can be represented in terms of an eigenvalue 
expansion. Recall that the potential qN is integrable on Ar. Let (Xk)keN be 
the sequence of eigenvalues of the operator A — qjy in Ar with zero boundary 
condition, and let (ek)keN be an orthonormal sequence of corresponding 
eigenf unctions. We assume that Ai is the principal eigenvalue and that e\ is 
positive in Ar. Then we have from (2.24) that 

E (N) [e -^( 9JV ,^ JV) > ] > e -|| /jv | U+{JV -l) inf v e (^+l)A* ( VNj e |) 

(2-25) keN 

> e -\\M\~+Mv e V +1 ^(u N ,et). 

Fix some bounded L 2 -normalized function G H 1 (M. d ) satisfying supp(/ijv) C 
Ar, then we may estimate 

Xi>-\\Vh N \\l-(q N ,h 2 N ) 

(2.26) >-\Nh N \\ 2 2 -(W-f N ,h 2 N )-\\h N \\lY.[ V$\dx) 

j<n JAr 

>-C r , 

where Cr is nonrandom and depends only on sup A R W, ||/at||oo) h]y and 
v (\x\) dx. Hence, we obtain from (2.25) that 

(227) EW[e- % ^>] 

1 ' ' > e -||/wll«.+(^-l)infw e -OK e -^[||Vh w ||g+(9 W ,h^>]/ I/Jw . )e 2\ 

The technical difficulty is now to find a positive lower bound for (z/jv, ef) that 
does not depend on the first N — 1 Brownian motions. Assume that R > 
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is so large that i^jv(-^il-i) > 0. We now use Harnack's inequality to obtain 
a pointwise lower bound for e\ in We first check that the random 

potential qw lies in the Kato class (see [8], Chapter 3, for more on the Kato 
class). It is clear that we only have to check this for the random potential 
V/i^p for j < N. We now verify that lim r |o VK r ) = 0> where 

(2.28) V(0= sup sup / G(x,y)V^{y)dy. 

neMi(R d ) xeA R J\y-x\<r 

This is seen as follows. First note that it suffices to show that 
lim sup / G(0, y)v(\y — w\)dy = 0. 

rlO weR dJ\y\ <r 

Recall that we are under the assumption J Be r ^ G(0, y)v(\y\) dy < oo, where G 
is the Green's function of the free Brownian motion, and v < v on (0,e), and 
v is decreasing. After possible alteration of e, we can extend v to a decreasing 
function (0, oo) — > (0, oo) such that v < v on (0,oo). We choose r G (0,e/3). 
We distinguish the cases \w\ > 2r and \w\ < 2r. In the first case, we have, for 
\y\ < r, that \y — w\ > r > \y\ and, therefore, v(\y — w\)< v(\y — w\) < v(\y\). 
Hence, 

sup / G(0,y)v(\y-w\)dy< [ G(0, y)v(\y\) dy, 

\w\>2r -> \y\<r J\y\<r 

and this vanishes as r { 0. 

In the case \w\ <2r, we split the integration over \y\ < r into the part 
where \y\ < \w\/2 and the part where |y| > |io|/2. On the first part, we 
have \y — w\ > \w\ — \y\ > \y\ and, therefore, v(\y — w\) < v(\y — w\) < v(\y\). 
Therefore, the first part can be estimated by 

sup/ G(0,y)v(\y-w\)dy< / G(0, y)v{\y\) dy, 

|«)|<2r^|y|<|w|/2 J\y\<r 

which vanishes as r J. 0. On the area where \y\ < r and \y\ > \w\/2, we can 
estimate \y — w\ < \y\ + \w\ < 3|y| and therefore find a constant C > such 
that G(0,y) < CG(0,y — w), and this means that we can estimate 

sup / G(0,y)v(\y -w\)dy 

\w\<2rJ\w\/2<\y\<T 

<C sup / G(Q, y — w)v(\y — w\) dy 

\w\<2rJ\w\/2<\y\<r 

<C [ G(0,y)v(\y\)dy, 

J\y\<3r 

and the proof of lim r jo ip(r) = is finished. In particular, this means that 
gjv is in the Kato class. 
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According to [8], Theorem 5.18, there is a constant A > 0, only depending 
on R and qw, such that 



(2.29) 



inf e\ > A sup e\. 
Afl-i Afi-i 



A closer inspection of the proofs of [8], Proposition 5.16, Theorems 5.17 
and 5.18, shows that the constant A does not really depend on the function 
qN, but only on the numbers sup^ fl W, ||/jv||oo and on the small-r behavior 
of the function if) defined in (2.28). Hence, A is nonrandom and does not 
depend on the first N — 1 motions. 

So we have (i/jy, ef ) > A 2 VN{h-R-i) sup Aj! _ 1 e 2 . Now we get that sup Afl _ 1 e 2 > 
(2|A^_i|)~ 1 . To see this, note that otherwise we would have j AC 
and, therefore, 



ef > 



(2.30) 



Ai> 

> 



-H/iv||oo + (iV- 
-||/iv||oc + (iV- 



l)miv + (W,el) 

1) inf v + i inf W — > oo as R 



oo, 



according to the assumption on W in (1.1). But, as we shall see at the end 
of the proof, this is impossible if is chosen appropriately. Summarizing, 
there is a constant C > 0, only depending on N, W, v, i>n and R, such 
that, almost surely with respect to the first N — 1 motions, 



E 



(N)\ 



131 W\ 



(2.31) 



\\Vh N \\l + (W- 



j<N 

We iterate now this argument for the ith motion for i = N—l,N — 2,...,l. 
For doing this, we have to replace the random potential gjy by qi = W — /j + 
J2j<i ^ 7*!g +Sj>i and obtain, for any bounded L 2 -normalized function 
hi £ H 1 (M. d ) satisfying supp(/ij) C A/j, almost surely with respect to the first 
i — 1 motions, 



(2.32) 



> Cexp<^ -P 



I 2 4- 



where C does not depend on (3 nor on the motions. This gives, recalling 
(1.27), 



liminf — logE[e 1 

P^OD (3 



N,0 



-K 



> 



N 

E 



N,0- 



N 



-</, 



v^ni + E< w - ^ + E^ ™?) 
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Now maximize the right-hand side over all choices of bounded L 2 -normalized 
functions h\,...,h^ £ H 1 (M. d ) satisfying supp(/i.;) C Ar. It is easy to see 
that, in the limit R — > oo, we obtain that the maximum of the right-hand side 
tends to —N\^\f)- Furthermore, one can see that it is possible to choose 
approximate maximizers hi,...,hN (depending on R) such that 
limsup^^^dl V/iilH + {W,hjy)] < oo and limsup^^ ||/ij||oo < oo. In partic- 
ular, the eigenvalue Ai introduced below (2.24) satisfies limsup^_ >00 (— Ai) < 
oo, and this explains why (2.30) is not possible [see (2.26)]. 

This completes the proof of the lower bound in (2.2) in the case of a 
soft-core potential v. 

2.4. Finish of the proof of Theorem 1.7. Proposition 2.1 implies the ex- 
istence of the logarithmic moment generating function A^y in (2.1) and 

identifies it with the Legendre transform of the function 1^ defined in 
(1.29). In particular, Theorem l.T(i) is implied. 

Now we prove the large deviations principle in Theorem 1.7(h). We use 
the Gartner-Ellis theorem (see [11], Corollary 4.5.27). For doing this, it 
suffices to show that the family of tuples (pp , ■ ■ ■ , ftp ) is exponentially 

tight under P^l as /? — ► oo, and that is Gateau-differentiable. The first 
condition is verified as follows. We follow the technique of the proof of [11], 
Lemma 6.2.6. Pick sequences Rk — > oo and [ and put K = H/cgm!/ 1 e 
Aii(M. d ) : /j,(A r ) < Sk}- The Portmanteau theorem implies that K is closed, 
and Prohorov's theorem implies that K is relatively compact, hence, K is 
compact. Now fix i S {1, . . . , N} and note that 




Now let a large R > be given. We additionally require that inf A flfc W > 
Rk for all k £ N [here we use our assumption in (1.1)]. We denote K now 
by Kr. Using the fact that W > and that v is bounded from below, it is 
easy to derive the existence of some constant CsK, not depending on R, 
such that limsup^^ J logP^(^ } g K R ) <-R + C for all R > 0, and this 
implies the exponential tightness. 

The Gateau-differentiability of is proven as follows. The proof of 
Lemma 1.3 shows that the infimum in the formula on the right-hand side 
of (2.2) is attained. We abbreviate A = A^f . Fix $ = / = {f\, . . .,f N ) £ 
Cb(R d ) N and some g £ C^(M. d ) N . We want to show the existence of the limit 
limt_»o t[A(/ + tg) — A(/)]. With (/4 , • • • i^n) a minimizer for the formula 
on the right-hand side of (2.2) for / replaced by / + tg, we obtain, by 
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replacing the minimizer in the formula for / by (/xf , . . . , pL^j ) , 

(2-33) -[A(f + tg)-A(f)]>J2(9i^?)- 

1 i=i 

Since the family (//^ , . . . ,/4v )t>o is easily seen to be convergent weakly 
toward the minimizer . . . ,[in) for the formula for A(/), it is clear that 
the right-hand side of (2.33) converges toward I^Li(ft; - Analogously, one 
shows the complementary bound. This implies the Gateau-differentiability 
of A { P with 

r) N 
(2-34) g-AjfV) = 

" i=l 

Now [11], Corollary 4.5.27, implies Theorem 1.7(h). The statement in 
Theorem 1.7(iii) is a standard corollary. This ends the proof of Theorem 1.7. 

3. Large deviations for the Dirac-interaction model: Proof of Theorem 1.12. 

We follow the same strategy as in the proof of Theorem 1.7 in Section 2. 
Hence, the main step is the proof of the following. 

Proposition 3.1 (Asymptotic for the cumulant generating function). 
For any f f N eC h (R d ), 

lim r p/{2+p) logEle-^L /„" 

/3— >oo 

(3- 1 ) x e -^* + * )/(a+rt Ei S1 < i < w 4 4,, V /C!H * )</ '' i ? > ] 

= -iY X 5J o) (A,/), 



iV#(A,/)= inf 



h 1 ,...,h Jv eH 1 (IR [i ): ||/ii||2 = lVi 

(3.2) 



J2h(hl) + (20-f,(h 2 ^ 

i=l 



l<i<j<N 



where we wrote (/i 2 )® = h{ © ■ ■ ■ © hr N and / = /i © • • • © /tv- 



Theorem 1.12 follows from Proposition 3.1 in the same way as Theo- 
rem 1.7 follows from Proposition 2.1 in Section 2.4; we omit the details. 
Hence, it remains to prove Proposition 3.1, which we do in the next two 



TRAPPED INTERACTING BROWNIAN MOTIONS 35 
sections. Let us first remark that 

W(sP)dt = f3 p / ( - 2+p \W,fif) and 

.(0 W..0"), 



(3 Jo 



f3 ( d+P )/(2+p) a (i,j) = pp/(2+p)/^P_^ d ^£_\ 

\ dx dx j 

as is derived by an elementary calculation. According to Lemma 1.16, ((J,g )/3>o 

satisfies large deviations principle bounds with scale f3 p /^ 2+p \ Hence, the 
problem left to be solved is to circumvent the missing boundedness and 
continuity of the functionals fit— > (W,n) and {ni,^) l— ► (~Jsr>~j5r)- 

3.1. Proof of the upper bound in (3.1). We proceed as in Section 2.1. 
However, an additional smoothing argument will be necessary. 

Let auxiliary parameters rj > 0, M > and i? > be given, recall Ar = 
[— R, R] d , and consider the A/j-periodizations Lg R of the densities lL de- 
fined in (1.37). As in Section 2.1, we distinguish whether or not ^(Ar) > 
1 — rj. Furthermore, we estimate 

(W,4>)>(W M ,4>) and (^,^)>(4<), i »AM>, 

where we recall that Wm = H^AMis the cut-off version of the trap potential. 

We intend to replace Ly by L\ R . The replacement errors for the terms 
involving W and / have been estimated in Section 2.1. The one for the 
interaction term is estimated as follows, 

(Lf,L^AM)-(Lf >R ,L%AM) 



> 



Lf{Lf A M){x) dx- L% R (x)(L% R A M)(x) dx 

[ ( L f ( x + 2 R k ) { L f A M) (x + 2Rk) 

k&Z d \{0} Ar 

- Lf tR (x)(L$ R A M){x)) d(x + 2Rk) 



>-M^f(A R )>-r,M. 
Hence, as in Section 2.1, we obtain the bound 

-/3»/( a -+*)(uinf A « :W-C f ) 
y > j_ e l3P^ 2 +P^ v [NM+N 2 M+Cf] 
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The last functional in the exponent on the right-hand side is bounded, but 
not continuous. Hence, we need a smoothing argument. For this purpose, let 
k$ denote the Gaussian density in M. d with variance 5 > 0, and denote convo- 
lution by *. Fix some small e > 0. According to [15], Lemma 3.7, specialized 
to our situation, 

(3.4) limlimsup/r p/(2+p) log sup P(| {if, g - g * k s ) | > e) = -co. 

5J.0 /j^oo ||s||oc<M 

Actually, in [15] only the discrete-time case is handled, but the continuous- 
time case is similar. This means that we can estimate the last expectation 
on the right-hand side of (3.3) from above against 

e -pp/( 2 +p)XC(M,s,S) 

where C(M,e,5) is a constant that satisfies lim^o C(M, e, 5) = oo for any 
M,e > 0. The functional 

^ ~ ( A M * **Htr «<"> d ~B AM * ^) 

is bounded and continuous in the weak topology on the set of probability 
densities on Ar. Hence, we can apply Varadhan's integral lemma and the 
large deviation principle in Lemma 1.16(i) to the expectation in (3.5). This 
gives that the limit superior on the left-hand side of (3.1) is not bigger than 



(3.6) 
where 



mmjr/inf W - C f , -rj[NM + N 2 M + C f ] + \C(M, e, S), 
- rjN[NM + N 2 M + C f ) - Xe - N Xn °\R, M, 5, A,/H, 



( N 

= inf \Y1 4 PCT) (to) + (Mm ~ f, Mi ® • • • ® »n) 

U=i 



. '. \ dx ' V dx 

Hi,..., hn G A^i(A^)|. 

Now we let 5 j 0, e [ 0, R — ► oo, 77 J. and finally M — ► 00. It is elementary 
to derive that limin%) Xn°\R, M, 5, A, /) > x N o) {R, M, 0, A, /), where we 
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interpret fi * kq as [i. Furthermore, similarly to the proof of (2.14), one 
shows that 

liminf liminf #(#, M, 0, A, /) > #(A, /). 

M— too H— >oo 

This ends the proof of the upper bound in (3.1). 

3.2. Proof of the lower bound in (3.1). For the proof of the lower bound, 
we follow a discrete variant of the strategy used in Section 2.3, that, we use 
iV eigenvalue expansions for the TV" expectations over the N random walks 
separately. First we consider the expectation with respect to the iVth walk, 
almost surely with respect to the other iV — 1 walks. We only treat the terms 
depending on the iVth motion and introduce the potential 



Qn(z) = £p' 



W 



(I 



j<N 



(I) 



Hi 



where = ^ 2+p \ and f N (x) = ^I x+[0>$ -y fs(y) dy. Denote by E« the 

expectation with respect to the ith random walk and by v\ its initial 
distribution. Choose R > so large that ^(-Bfii-i)^) > 0. Then we have 



E 



(3.7) 



> 



£ u N (z)Ef)[e-Io ^^ (W) )*l{ B upp(4 JV >) C B^jtW 



(N) 



> VN(z)e N (z) 2 e 



I3\n 



z£B 



where A at is the principal eigenvalue of A — in Br£„ with zero boundary 
condition, and e/v is the corresponding positive ^-normalized eigenf unction. 
Pick some L 2 -normalized function hjy € C 2 (K d ) satisfying supp(ZiAr) C Ar 
and pick vn(z) = (1 + o(1))^ /3 ^ 2 /iat(z^ 1 ), where o(l) is chosen such that 
vn is £ 2 -normalized. Then we have, using the Rayleigh-Ritz principle for 
the principal eigenvalue, and noting that f N — > /jv as (5 — > oo uniformly, 



Atv> 



\ j<N I 



X (1 + 0(1)) 



II VMI 2 + {W- /at, h%) + A £ (Lfihl) 

j<N 
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We estimate the term J2z v n{z)zn{z) 2 on the right-hand side of (3.7). Note 
that eTv is also an eigenfunction for the transition densities of the random 
walk in -Br^ with potential — q^v -Aat, that is, 

e N (z) =E*[e-/oW*)+AAr)* 1{5t G Vt G [0, l]}e 7V (5 1 )], 

where E z is the expectation with respect to an independent copy (St)te[o,oo) 
of, say, (S\ ; )te[o,oo)- We can estimate Atv < £1 2 ||/tv ||oo and, since Ly(x) < 
£j for any x G M d and any j G {1, . . . , A/"}, 



9iv(5 t ) < £ 







R P + \Y. L f( S t) + WfN\\oo 

j<N 



< [R P + ||/jv||oo]/r 2/(2+P) + X(N - 1)/^-2)/(2+ P ) < C(3 (d-2)/(2+p)^ 

for some C > 0, depending only on i?, A, iV, p and ||/jv||oo- Using this in 
(3.8), we find, for all large (5 (changing the value of C if necessary), 

e N (z)>e- cl3(d - 2)/i2+P) P^teB^Vte [0,11,5-1 = 2)6^(5), 



z G Bp^p . 

Recall that p > d - 2, hence, /3( d " 2 )/( 2 +P) = (/3?V( 2 +p)). It is clear that 

P 2 (5 t G 5 flC/3 V< G [0, = z) > e -°(/3 p/(2+p) ) 5 uniformly in z, z G Bjj^. Since 

e^v is ^-normalized, we can estimate J2j^B Ri e A r ( z ) — E^gs^ 

Pick some z^v G -B^ such that un(zn) > 0, then we have from the preced- 
ing, for any large f3 > 0, 

E W[ e -V/3/ V(5 t W )^-A/3(^)/(2+P)^ j< ^ a ^+ / 3P/(2+p)( / ^,LW) ] 

> expl -/?p/( 2 +p) 



(3.9) 



IIV^Hl + - /at, 4) + E <^X) 

+ ( /3 p/(2+p)) 



and the term o(/? p ^ 2+p -*) does not depend on the iV — 1 other random walks. 

In the same way, we treat the expectations with respect to the other 
N —1 random walks. For doing this, introduce recursively, for i = N —1,N — 
2, . . . , 1, the potentials 



*0) =ip 



j<i 



j>i 
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where fi(x) = £^ / x+ [o ^ 1 ) d h(v) dy, and pick L 2 -normalized functions Tijv-Ij 

hjsi-2, ■ ■ ■ ,hi G C 2 (M. d ) with supports in Ar. In the same way as (3.9), one 
derives 

Ef e -V/'Eti/>(^)V^ ,/M ^<*<i^ ^V /1 ^ ) ^ > l 



>e°^ /(2+p) )exp J-/5p/(2+p) 



A? 



E0l v ^ll2 + (^-/iA 2 >] 



-a E to a >*3 

l<i<j<N 

Picking the L 2 -normalized functions h\,...,htf with support in Ar opti- 
mally and letting R — > oo, we arrive the lower bound in (3.1), noting the 
substitution hf(x)dx = fii(dx) in (3.2). 

4. Analysis of the ground product states. We prove Lemma 1.3 in Sec- 
tions 4.1-4.4 and Theorem 1.14 in Section 4.5. The proofs combine methods 
from variational analysis and from probabilistic potential theory. 

4.1. Existence of minimizers in (1.11). Let (hi^, . . . , /ijv,fc)feeN be a se- 
quence of approximate minimizers for the right-hand side of (1.11), that is, 

hi t k € H 1 (M. d ) with \\hi t kh = 1 for all t = 1, . . . ,-ZV and Jfe G N, and Nx^ = 
lim/^oo^^Li || V/ij^Hl + (2H+ 0, hfy], where we abbreviated h k = hi t k®- ■ -® 
tiN,k- Fix i G {1, . . . ,N}. Our first goal is to establish the convergence of 
(^i,fc)fceN along suitable subsequences. 
For any R > and any fcEN, we have 

9 , s Sin* h? h (x)W{x)dx 

hl k {x)dx< m l : k ) ' ) 



R 



(4-1) < 

< 



inf A ^ W 

E£i II VMIl + (2H + 0, h 2 k ) - ((N(N - l))/2) inf t, 



inf A c fl 

By our assumption in (1.1), sup fcgN J A c hf k {x) dx vanishes as R — ► oo. Hence, 



the family of probability measures {hf k {x)dx) k ^ is tight, and from Pro- 
horov's theorem, we conclude the existence of a probability measure //j G 
.Mi(R d ) such that (hf k (x) dx) ke ^ converges weakly toward /ij along a sub- 
sequence, which we again denote (k)keN- Hence, for any bounded continuous 
function / : ~R d — ► R, we have lim/ c ^ 00 (/, /i? fc ) = (/, /ij) . 

Furthermore, because H^/d^ = 1 for any k G N, according to Banach- 
Alaoglu's theorem, there is a hi G L 2 (R d ) such that {hi :k ) ke ^ converges 
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weakly in L 2 along a subsequence [which we again denote (fc)fc e pj] toward hi. 
Since, by approximative minimality and because W and v are bounded from 
below, also the sequence of energies (||V/ti ) fe||2)fegN is bounded, the Banach- 
Alaoglu theorem implies that also the sequence of gradients converge weakly 
in L 2 along a subsequence [which we again denote (fc)/. g pg] toward some gi. In 
particular, hi £ and V/ij = gi, which follows via [21], Theorem 8.6, 

from the completeness of the Sobolov space. By [21], Theorem 8.7, we may 
assume that the convergence h^k — > hi is pointwise almost everywhere. 

We argue now that h 2 is a density of //j. According to [21], Theorem 8.6, 
the convergence of (hi^)k toward hi is also strong in L 2 on every compact 
set. In particular, the integrals of h 2 k against any continuous function with 
compact support converges to that of h 2 . By weak convergence of h 2 k (x)dx 
toward /ij, they converge also toward the integrals against /ij. Hence, the in- 
tegrals of any bounded, compactly supported function against hf and against 
fj,i coincide, which implies that /ii(dx) = h 2 (x) dx. 

Now we prove that the vector (pi, . . . ,[in) is a minimizer on the right- 
hand side of (1.11). To do this, we verify suitable lower-semicontinuities of 
the functional on the right-hand side of (1.11). It is well known that the 
energy functional, h>— > ||V/i||2, is weakly lower-semicontinuous in 1? . Fur- 
thermore, the map fi (W, /u) is obviously lower-semicontinuous in the weak 
topology of probability measures, since W is nonnegative and continuous. 
By c + and D_ we denote the positive and the negative part of d, respectively. 
Fatou's lemma implies that liminffc^ 0O (d+, h\) > (d+,fo 2 ), where we recall 
that hk = h\ t k ® ■ • ■ ® hjy^- Furthermore, note that 

\(* + ,h 2 )-(» +1 h 2 )\ 




Using (4.1) and the strong L 2 convergence on compacts, one easily derives 
that the right-hand side vanishes as k — > oo. Summarizing, 

N 

V||V^||2 + (2B + D,/i 2 ) <liminf 

k— »oo 

i=l 

where we wrote h = h\® ■ ■ ■ ®h^ and (in an abuse of notation) hk = h\k <g> 
• • • (£> h^,k- Hence, the tuple (hi, . . . , /ijv) is a minimizer of the right-hand 
side of (Lll). 

4.2. Positivity and regularity of minimizer s, soft-core case. First we con- 
sider the case of a soft-core potential v. Let (hi, . . . , /ijv) £ i/ 1 (M a! ) Ar be a 
minimizer of the variational formula in (1.9). Fix i £ {1, . . . ,N}. We prove 



/v 



J2\\vh k \\ 2 + (<m+»,h 2 k ) 



N 



i=l 
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that hi is strictly positive everywhere in M. d . We certainly may assume that 
hi > 0. Note that hi is a minimizer in the problem 



(4.2) Xi 



inf 



|^i||2 = l 



3IIVM2 + <w, /*?> + ( hiY^vh) 



First we show that Vhj is locally integrable for any j G {1, . . . , N} with j ^ i. 
Indeed, for any measurable set B C M d of finite measure, we have [denoting 
the e-ball around x by B £ {x)] 



Vhj(x)dx = J dyv(\y\) J dxhj(y + x) 



(4.3) 



< 



< 



B e (0) 



[e,oo) 



dyv{\y\)\\hj\\l + sup v / cte / dyh 2 (x + y) 



B 



B|(0) 



^y^(|y|) + l-^l su p ^ 

Be(0) [£,00) 

< 00. 



Since Vhj is also bounded from below, the interaction potential of the prob- 
lem in (4.2) is locally integrable. The same holds for the trap part, thus, the 
potential W + Y^j^iVhj ~ i n (4-2) is locally integrable. By [21], Theorem 
11.8, hi satisfies the Euler-Lagrange equation 



(4.4) Ahi = Uihi in V'( 



where Ui = W + J2 vh ) ~ x i- 



(Note here that the potential Ui is a priori not locally bounded, but the first 
part of the proof of [21], Theorem 11.8, does not use this.) So, we get by 
[21], Theorem 11.7, that 



(4.5) 



for all q < 



00, 
3, 



in d = 2, 
in d = 3. 



Now proceed for d = 2 and d = 3 separately. In d = 2 we deduce that 
Vhj G L^ oc (]R 2 ) for any p G [1, 00). Indeed, fix any bounded measurable set 
B C M. d , pick rj > so small that v(\y\) > for all y G B v (0) and abbreviate 
Krj = f B , ->dyv(\y\). Then, using Jensen's inequality and (4.5), we see that 



B 



Vh 2 Ax)\ v dx<K v \ dx 



(4.6) 



B 



< K J B dx 



+ \B\ sup \v\ p 
(77,00) 



B v (x) "77 

d y V ^ X - y ^ h j {yfP+\B\ sup \v\P 



B v (x) 



(77,00) 

RP- 1 ( dyv(\y\) f dx h 2p (y + x) + \B\ sup \v\ p 
Jb v {o) Jb (77,00) 



< 00. 
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Hence, the potential U is locally p-integrable for any p < oo. Using again 
[21], Theorem 11.7, we obtain that hi is continuously differentiable, and 
all first partial derivatives of hi are Q-H61der continuous for any a < 1. In 
particular, C/j is also locally bounded, and from [21], Theorem 9.10, we have 
hi > everywhere in R 2 . For this, we apply [21], Theorem 9.10, to any open, 
bounded and connected j4cK 2 , that satisfies J A hi(x) dx > 0. 
Now we turn to the case d = 3. Recall that we assume that 



Bi(0) 



\v (\y\)\ 1+s dy < co for some 5 > 0. 



We show that |U/i 2 | p is locally integrable for some p > 
J~hf'(x)dx< 



for any p' < |. Pick a bounded open ball -B C 



Recall that 
^ 3 such that 



)gh^ p (x) dx < co. Pick some open ball -B whose closure is contained in B. 
We shall show that J B |U/i 2 (x)| p dx < co for some p > §• Let now 77 > be so 
small that B v (x) C B for any x £ B and that u(|y|) > for all y £ B v (0). We 
estimate |V/i 2 (x)| < sup^ >00 ) \v\ + Jb (o) dy v{\y\)h?-(x + y) for x E 5. Hence, 
we have to show that 



(4.7) 



dx 



B 



B„(0) 



*M|y|)/i 2 (x + y) 



< oo for some p > 



3 
2 ■ 



We shall do that for any p > 3/2 and p' < 3/2 satisfying p/p' < 1 + 5. Recall 
the abbreviation K v = j B , Q s dyv{\y\). Jensen's inequality gives 



dx 



B 



B„(0) 



dyv{\y\)h 2 Ax + y) 



< K p 



dx 



B 



K p-p/p' 
J1 tj 



B„(0) 
c/.r 



x + y) 



p/p' 



fcf (x + y) 



(4.8) 



B„(x) 



,p/p' 

cfe/i 2p (z) 



«(|y|) 



p/p' 



< TSP-P/P' 
— ±i -r) 



dx 



B 



B n (x) 



B„(0) 



dyfcf (x + yMlyl)^ 



.7 

P/P' 



< OO, 



where we used that p < p'(l + 5) in the last step. Hence, Ui 6 Lf 



Ioc 



for 



some p > 4. Now [21], Theorem 11.7, implies that hi is continuous and 
that Vh? is locally bounded. Therefore, Ui is locally bounded. Finally, we 
get from [21], Theorem 11.7, that hi is continuously differentiable and that 
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the partial derivatives are a-H61der continuous for any a < 1; whereas [21], 
Theorem 9.10, gives that hi is positive everywhere in R 3 . 

4.3. Positivity and regularity of minimizers, hard-core case. Now let v 
be a hard-core potential, and recall the parameter a G [0, oo) from (1.2). Let 
(hi, ■ ■ ■ ,hw) be a minimizer of the formula in (1.11). Fix i G {1, . . . , N} and 
put 



£li = B a [ |Jsupp(/i j ) 

Wi 

where B a (A) is the a-neighborhood of a set A if a > 0, and Bq(A) = A. Then 
it is easy to see that, for any j / i, the potential Vh 2 is bounded on any 
compact subset of fij. It is clear that supp(/ij) C fij. Hence, /ij is a minimizer 
in the problem 

(4.9) _ inf 

fe l e-ff 1 (K d ),supp(/i i )cO i ,||h i || 2 =l 

Again by [21], Theorem 11.8, hi satisfies the Euler-Lagrange equation in 
(4.4) in T)'(Qi). By [21], Theorem 11.7, hi is continously differentiable in Qj, 
and its first partial derivatives are a-H61der continuous for any a < 1. By 
[21], Theorem 9.10, hi is strictly positive throughout any open component 
of Qi in which it is not identically equal to zero. 

4.4. Boundedness of minimizers. Now we prove that ||ftt||oo ^ Cd(Xi — 
(N — 1) inf v) dj/4: . We do this by carefully examining the proof of [3], Corol- 
lary 2.5. Recall the potential 17% = W + ^j^iVh 2 j — A, from (4.4), which is 
locally bounded and globally bounded from below, more precisely, — inf 17 i < 
Xi-(N- 1) inf v =: m G [0, oo). Fix x G R d and < r < (4Cu i )- 1 / 2 , where 
C = sup| a .| <1 Ea;(Ti(0)), and T r (y) denotes the first exit time of a Brownian 
motion from the ball B r (y). Then we have ¥. XQ (2uiT r (xo)) = 2uiCr 2 < ^ and 
Khas'minskii's lemma ([3], Theorem 1.2) implies that 

(4.10) E^expj-^ X ° 2Ui(B s )d s yj < E Xo (e 2uiT ^) < 2. 

Note that hi(x) = E x .[e"/o Tr( " 0) u ^ B{s))ds hi(B(T r (x )))] for x G B r (x ) by [3], 
Theorem A. 4.1. We apply the Cauchy-Schwarz inequality and (4.10) to get 
that 

(4.11) h l (x ) 2 <2E X0 [h 2 (B(T r (x )))]=2 f a(dx)h l (x) 2 , 

J\x-xo\=r 




vmI + (w,h 2 ) + UlY^vh) 
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where a denotes the normalized surface measure on dB r (xo). Averaging the 
left- and the right-hand side of (4.11) on r € [0, (4Cuj) _1//2 ], we obtain 

hi(x ) 2 <2(ACui) d / 2 f hi{xf dx < 2(AC[X l - (N -l)mfv]) d/2 , 



and this implies the assertion. 

This completes the proof of Lemma 1.3. 

4.5. Proof of Theorem 1.14. Recall that we are in d S {2,3} and assume 
that v > and v(0) > and a(v) < oo, where a(v) is given in (1.41). Fur- 
thermore, recall that we replace v in (1.9) by vn(-) = N d " 1 v(- N). 

First we show the upper bound, limsup^r^go — Let <fi denote 

the minimizer in the Gross-Pitaevskii formula (1-40) with a = a(v). Picking 
(h\, . . . , /ijv) = (</>,...,(/>) in (1.11), we obtain that 

+ i / N d v{N\x\)((j) 2 *cp 2 )(x)dx 

JWL d 



(4.12) 



t + (w,4> 2 ) 



\ lv{\x\){cl> 2 *<t> 2 ){x/N)dx, 



where h * g(x) = J Rd h(y)g(x + y) dy. Since (j) 2 is bounded, integrable and con- 
tinuous, we can use Lebesgue's theorem to see that the last term converges 
toward \f v{\x\) dx{<j) 2 *^ 2 )(0) = 4vr5(z;)||(/)|||. Therefore, we have 



limsupxS? < ||V^|| 2 + (W,4> 2 ) +4ira(v)U\\j = X ^ . 



Now we show the reversed inequality, lim inf jv— ►oo Xn — X~/ v \ ■ Let, for 
any JV£N, (^-i^, • ■ • > hffl) be any minimizer in (1.11), and define 

N 



^4E(^ (a °) 2 - 

i=l 



Since, for every R > 0, 

(4 ,3) t , %(x)dx <^M<^ 



mi A c R W ~ inf A ^ W 

the sequence of probability measures (4> 2 N (x) dx) n& f>] is tight [recall (1.1)]. 
Hence, there is a limiting probability measure of 4> 2 N (x)dx in the weak 
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sense, along some subsequence. Since (4>n)n is bounded in L 2 (R. d ), there is 
an L 2 (M d )-weakly converging subsequence toward some <f) G H (R ). Since 
also (||V^>Ar||2)7v is bounded [which is derived similarly as in (4.13) with the 
help of the convexity of the map (p 2 t— > || V<^||2], we may assume that VcpN 

weakly in L 2 (M. d ). Furthermore, we may assume that the convergence 
(pN — > <P is even strong on every compact subset of M. d and pointwise almost 
everywhere. In particular, all integrals of <pM against continuous, compactly 
supported functions converge to their integral against 4>. Hence, ^(dx) = 
4>(x) 2 dx, and we have that \\(j)\\2 = 1- By L 2 -weak lower semicontinuity and 
convexity of the energy, we have 



1 N 

HV^Hl < Hminf IIV^HI < Uminf — ^11 V/i. (JV) l|2 

N— >oo N—too iv — * 

1=1 



and 



(W,(j> 2 ) < liminf (W,<&>. 



Hence, we have 

liminf ±-(h[ N) ® • • ■ ® h^\n N h[ N) ® • • • ® h { p) 
N—>oo ly 

> ihminf f «(M)($r?$r)(-|) dx+ \\V^\\ 2 + {W^ 2 ) 

-Ihmsup-LW f hf l \y) 2 hf l) (y+^)v{\x\)dxdy. 
Hence, for proving that lim inf n-^oo Xn — X~^> h remains to show that 



liminf / v(\x\)(<f> 2 N *<i> 2 N )(^-)dx 



(4.14) > / v(\x\)dx\ 

JM. d 

1 N f f ( X \ 2 

limsup— _ l, h i N) (y) 2h i N) \y + M) v (\ x \) dxd v 

(4.15) =0. 

For the proof of (4.15), recall from Lemma 1.3 that H^IU < C d (\[ N) ) d / 4 , 
where Cd > depends only on d, and A^ is the Aj of Lemma 1.3(h) with 
v replaced by vn- Hence, 

i N r r , v 2 
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1 N r 



i=l 



Recall that jj £^ 



xS?> + ^<M N } ) 2 ' M^)? ^ which 



is bounded from above in N, as we have seen in the first part of the 
proof. Hence, (4.15) directly follows in d = 2, and in d = 3 we estimate 

(A 4 (A °) 1/2 < CN l l 2 for some C > and all iV G N, and (4.15) also follows. 
We turn to the proof of (4.14). Fatou's lemma gives that 



lim inf 

N— »oo 



> 



,«(M)WW*) 



da;«(|x|) / dylimmf4> 2 N (y + ^-)(j) 2 N (y). 



Recall that limjv^oo 4>n( x ) = (j>(x) for almost every x G 
suffices to show that 



Therefore, it 



(4.16) lim sup 

]V->oa 



<P 2 N(y)-^ N {y + ^j 



for every i,t/£ 



Fix x,y G For a while, we write z instead of jr. Estimate 



1 



A? 



(y)-^( y + z )\<±\2\€ ] (y)-h { f , (y + z)\2\\h l 



(4.17) 



i=l 



N 



<^^rf A \hf\y)-hf\y + z )\. 



i=l 



In the following we denote by C a generic positive constant, which may 
change its value from appearance to appearance and does not depend on i 
nor on N. Fix i G {1, . . . , N}. By U\ ■ , we denote the potential in (4.4) with 
u replaced by vn, that is, 

U^ N \w) = W{w)+Y J V N h 2 iw) - \f\ 

where Vn is the operator V with v replaced by vn- Note that 

1 N r 
<T7E SU P djM>(JV|s - y^y) 



(4.18) 



< 



a 



2 N 



N 



V(Af } ) d/2 / V (|z|)^<CiV( d - 2 )/ 2 



Recall (4.2) and estimate 
(4.19) 



inf [||V^[|| + <W, /if) + CTV^- 2 )/ 2 ] 

h i eH 1 Qs. d ),\\hi\\ 2 =i 
< CAf (d ~ 2)/2 . 
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We write T :=T\(y) for the exit time of a Brownian motion (Bt)t>o from 
the ball B := B\{y). By [3], Theorem A. 4.1, Lemma A. 4. 5, we have 



hf ) (w)=E w [hf ) (B T )] 
(4.20) r ft 



/ UrHB s )hf\B s )ds 
Jo 



w€B. 



The first term is harmonic. By the mean- value property, we see that, denot- 
ing by S y>z the symmetric difference B 1 u(y)ABiu(y + z), 

\E y [hf \B T )] - E y+z [hf \B T )]\ < — 1— f \E w [hf\B T )]\dw 

\Bi/i{y)\ Js y , z 

(4.21) Hoc 

<C{>^) d ' A \z\. 

The second term on the right-hand side of (4.20) may be written as hff^ (w) + 
h[ N \w), where 



(4.22) hf\w):= G(w,w)U- N \w)h\ N \w) dw 

Jb 

(4.23) hf\w):= [ (Gb — G)(w, w)U- N \w)h\" 1 (w) dw, 



IB 

and G and Gb are the Green's function in M. d and in the ball B, respec- 
tively. (For explicit formulas for Gb, see [8], Section 2.3, e.g.) By [21], The- 
orem 10.2(h), for p = d and any a 6 (0, 1), we have 

By (4.18) and (4.19), sup weB \U^ N \w)\ < CN^ d ~ 2 ^ 2 . Hence, we can further 
estimate 

(4.24) \hf\y) - hf\y - z)\ < C\z\<* N^W 2 (xfY* ■ 

Now we turn to h\ N \ which is harmonic in B, since (Gb — G)(-,w) is 
for any w e B. By [8], Theorem 2.5, (G - G B )(w,w) = E w [G(B T ,w)] for 
any w,w £ B. Like in (4.21), we use the mean-value property, recall that 
su Ptugs \ Ui N \w)\ < CN( d ~ 2 ^ 2 , use Harnack's inequality and obtain, for any 
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ae(0,l), 

\hrHy)-hf\y-z)\ 



< CN (d-2)/2^N ))d/4 >_ / dw / M{G _ Gb){w ^ 

\ B i/i{y)\ Js y , z Jb 



(4.25) KCN^I'^f 4 — * / dw / dwE [G(B T ,w)] 

\tfi/4{y)\ Js y , z Jb 

< CN^WQtWf* , j%L , / G(0,«0dti5 

l^i/4 (y) I Vb 4 (o) 

Substituting the estimates (4.21), (4.24) and (4.25) in (4.17), we obtain 

i ^ 

my) - + *)i < c± ^(xfY 2 ^ + iv (d - 2)/2 i.r + N^\zn. 

i=l 

Recall that z = x/N and, furthermore, that jfJ2iLi^\ N ^ < C and < 
CN^ d ~ 2 ^ 2 . Picking a sufficiently close to 1, (4.14) easily follows. This fin- 
ishes the proof of limjy^oo = 

The proof of the convergence of </)n toward the Gross-Pitaevskii mini- 
mizer <j) (not only along subsequences) is done in a standard way as follows. 
Let / :M. d ->t be a measurable bounded function, and fix some e£R with 
|e| small. Consider the trap potential W e = W + ef. It satisfies all the re- 
quired assumptions that we posed on W, with the possible exception of the 
nonnegativity. However, an examination of our above proofs shows that they 
work also for W e in place of W. Let <Mr(e) denote the variational formula 
in (1.11) with W replaced by W e , then we know from the preceding that 
limAr^oo gAr(e) is equal to the Gross-Pitaevskii formula in (1.40) with W 
replaced by W e . Furthermore, it is easily seen that is concave in a neigh- 
borhood of zero for any N. Hence, we know that also the derivatives g' N 
converge toward the e-derivative of the Gross-Pitaevskii formula with W 
replaced by W e . Clearly, g' N (0) = {f,4> 2 N ), where (f)N is as in the above proof 

of liminf jv->oo Xtv*^ — ^Sy) Furthermore, the analogous derivative for the 

Gross-Pitaevskii formula is equal to (/, cf) 2 ), where <f> = 0—/^ is the min- 

imizer in (1.40) for a = 5(f). By convergence of the derivatives, we have 
liniAr_ >00 (/, 4> 2 N ) = (/, (j) 2 ). Indeed, for e > 0, we have 

limsup(/,(/) 2 v ) < limsup-(c/Ar(e) - ^(0)) = -{g(e) - g(0)), 

and the right-hand side converges to (/, 4> 2 ) as e [ 0. In order to see the 
reversed inequality, replace / by — / and e by — e. This shows the conver- 
gence of 4>%{x) toward (j)(x) 2 both in the weak L 1 -sense and in the sense of 
probability measures. 
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5. Large deviations for the canonical ensemble model: Proof of Theo- 
rem 1.5. In this section we prove large deviations statement on the canon- 
ical ensemble model in Theorem 1.5. We follow the same strategy as in the 
proof of Theorem 1.7 in Section 2. Hence, the main step is the proof of the 
following. 

Proposition 5.1 (Convergence of the logarithmic moment generating 
function). For any f £ C\,(Q), 

(5.1) lim ]-logE[e- l3 ^~- f ^] = -N X N(f), 

p— »oo p 

where 

(5.2) XN (f) = l inf [|| V / l ||2 + (2IJ + D -/,/ l 2 )]. 

iv heH 1 ^) : \\h\\ 2 =l 

Proof. See Sections 5.1 and 5.2 for the upper and lower bound, respec- 
tively. □ 

From Proposition 5.1, one derives Theorem 1.5 in the same way as The- 
orem 1.7 is derived from Proposition 2.1 in Section 2.4; we omit the details. 

The assertion of Proposition 5.1 is classical and well known if the potential 
2B + V — f were bounded and continuous. But 2H and d are unbounded to 
oo, and this is the additional technicality we are facing. Nevertheless, we feel 
that this problem and its solution are well known, but we could not find a 
precise reference. Hence, we provide a proof. 

We shall employ eigenvalue expansions on large compact sub boxes to 
derive (5.1). For the upper bound, we use periodic boundary condition in 
the box, and for the lower bound, we use zero boundary condition. Let 
us remark here that, if W and v would have been assumed continuous in 
{W < oo}, respectively in {v < oo}, these arguments for the upper bound 
could have been also replaced by applications of Varadhan's lemma and the 
large deviations principles in Lemma 1.15. However, due to the degeneracy 
of v at zero, we did not find any way to derive the lower bound in (5.1) via 
large deviations arguments in the soft-core case. 

5.1. Proof of the upper bound in (5.1). Recall that A R = [-R,R] d . We 
shall divide the path space into the part on which the Brownian motion 
spends more than (1 — r))/3 time units in up to time j3 and the remaining 
part where this is not satisfied. We will replace up by its A^-periodized 
version and control the error. We also replace the two potentials W 
and v by their cut-off versions Wm = W AM and % = v A M, where M > 
is large. Finally, we let R — ► oo, rj J, and M — ► oo. 
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In the following we sketch only the proof because all the details can be 
found in the corresponding proof for the Hartree model. For the canonical 
ensemble model, it turns out the error estimates are simpler, in particular, 
the interaction part, because here only one Brownian motion is involved, 
and hence, there is only one time scale. We begin with 



(5.3) 



E( e -W-/.^>) < ECe-^^lO^CAg) < 1 - r}})e- pMv eP^ 
+ M( e -^ WM+t ' M - f '^h{fJi l 3(A%) > 1 -r]}), 



where 2Um and X>m are defined as 2U and t> in (1.3) with W and v replaced 
by W A M and v A M, respectively. The first term on the right-hand side 
of (5.3) is easily further estimated, using that W > 0. Indeed, it is not bigger 
than exp{— /3(?/infA^ W + inf v — ||/||oo)}- For estimating the second term, 
we replace fip by its periodized version fJ>/3,R as in (1.45) with the box A 
replaced by A^. Now we can estimate the replacement error for the trap 
part as follows: 

N 

= E / jv E( w (^+ 2 ^)-^)W(^+2^)) 

^ b > keZ dN \{0} A R i=l 

> -M^{{K N R f) 

> —TjM. 

The replacement errors of the other parts of the potential are estimated in 
the same way: 

(t>Af - f,Hp) ~ {*>M - f,V/3,R) 

>(-M + inf V -2||/|| 0O V /3 ((A^) c ) 
= -r ) {M-miv + 2\\f\\ 00 ). 
Summarizing, we obtain 

l.h.s. of (5.1) < e - PivMA n W + M *-\WM 
{ 1 + e 2f3 V (M-ird v+\\f\\ ao ) E / e -P(Wu+*M-fM,ii)\ 

Now we use an eigenvalue expansion to derive that 

(5.6) limsupilogE( e - /3(2ITA/+t,A/ - / ' /i ' 3 - R) ) < -N X n(R,MJ), 

P^oo P 

where 

X n(R,MJ) = ±-m{1 f \V R h(x)\ 2 dx + (W M + *M-f,h 2 ): 

(5.7) N ^Agnn 

hec 00 (A%nn),\\h\\ 2 = i\. 
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Here Vr denotes the periodized gradient on A|j, that is, the one with peri- 
odic boundary condition. To derive (5.6), we let (Xk)ke'N an d (e^ken be the 
sequence of eigenvalues and an orthonormal basis of corresponding eigen- 
functions of the operator A — 2$^,/ — x>m + / hi n £1 with periodic bound- 
ary condition. We may assume that Ai = —NxN(R,M,f) is the principal 
eigenvalue. Let p\ (x,y) denote the transition probability function of the 

Brownian motion on the torus n O and note that p^\x,y) < Cr for 
any x, y E A^, for some Cr < oo. For simplicity, we assume that i^(A^) = 1. 
Then we can estimate 



= e~ inf f dy f u{dx)pf ] (x, y)E y (e-^- 1 ^ WM+ ° M ~^- 1 ^ ) 

JA^ J 

<Cj2e^~ 1)X Hek,t) 2 , 

fceN 

where C > does not depend on /?, and (■, •) is the inner product on L 2 (A|j). 
Now use Parseval's identity to continue 

E(e -/3<2H M +»M-/,^,fl>) < gOS-lJAx J2(e k ,l} 2 

= e 03-l)Ai|| 1 ||2 j 

where || • || denotes the I? norm in A^. From this, (5.6) directly follows. 
From (5.5) and (5.6), we obtain 

l.h.s. of (5.1) 

( 58 ) < -min|?7inf VF + infu - H/lloo, 

- 2 V [M -m£v+ H/lloo] + N XN (R, M, /) 

Now we let R — > oo, then rj j and finally M — > oo on the right-hand side 
of (5.8). Note that 

(5.9) liminfnminf X jv(i?,M,/) > X nU')- 

M— »oo R—>oo 

The proof of (5.9) is standard, we do not carry it out here; note that the 
proof of (2.14) is similar and even technically more difficult. 

Because of our assumption on W in (1.1), lim^^oo inf a c W = oo. Hence, 
we obtain that the left-hand side of (5.1) is not bigger than Nxn(I)- This 
completes the proof of the upper bound in (5.1). 
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5.2. Proof of the lower bound in (5.1). Now we prove the lower bound 
in (5.1). Like in the proof of the upper bound, we rely on an eigenvalue 
expansion. We recall Ar = [— R, R] d and the set U v from (1.5). Fix some 
rj > a [recall (1.2)]. We estimate the expectation on the left-hand side of 
(5.1) by imposing zero boundary condition in a certain compact set: 

E [ e -0<ar+i>-/,A*>] > E[e-PV° +0 -f>wh{supv{np) cA^n U v }], 



where we put U T) = U„ in the hard-core case and U n = K in the soft-core 
case. Hence, in the hard-core case, we may replace 2B and D by TSm and 
Vm [for this notation, see below (5.3)], respectively, where M > depends 
only on r\ and R. In the soft-core case, the potential 2U + V — f is locally 
integrable in the box A# and bounded from below. 

Consider the linear operator A — 2U — d/ on the space L 2 (A^ nU v ) with 
zero boundary condition. Standard results imply that this operator possesses 
a compact resolvent. Hence, a Fourier expansion in terms of the eigenvalues 
of this operator yields that 

(5.10) ^logE[e-^ +D -^)l{supp(^) C A%nU v }] 

where 

xn(R,vJ) 

(5 - U) =-1; mf _ (hvh\\ 2 2 + (W + V -f,h 2 

iv heH 1 (m dN ),su P p(h)cA%nU v ,\\h\\ 2 =l \ z 

is the principal eigenvalue of the above operator. The proof of (5.10) is 
similar to the one of (2.31) and technically much easier; we omit the details. 
Now we let R — > oo and n [ a and see easily that 

lim HmxN(R,vJ)=XN(f)- 

H— >oo r/la 

This completes the proof of the lower bound in (5.1). 
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